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We consider, from several complementary perspectives, the physics of confinement 
and deconfinement in the 2+1 dimensional Georgi-Glashow model. Polyakov's 
monopole plasma and 't Hooft's vortex condensation are discussed first. We then 
discuss the physics of confining strings at zero temperature. We review the Hamil- 
tonian variational approach and show how the linear confining potential arises in 
this framework. The second part of this review is devoted to study of the de- 
confining phase transition. We show that the mechanism of the transition is the 
restoration of 't Hooft's magnetic symmetry in the deconfined phase. The heavy 
charged W bosons play a crucial role in the dynamics of the transition, and we 
discuss the interplay between the charged W plasma and the binding of monopoles 
at high temperature. Finally we discuss the phase transition from the point of 
view of confining strings. We show that from this point of view the transition is 
not driven by the Hagedorn mechanism (proliferation of arbitrarily long strings), 
but rather by the "disintegration" of the string due to the proliferation of branes. 
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1 Introduction 



Study of the confining regime in QCD is one of the most elusive and compli- 
cated subjects in modern particle physics. The theory is strongly interacting, 
and our arsenal as of today does not contain any reliable analytical tools to 
tackle strongly interacting quantum field theory. The situation is exacerbated 
in QCD by the fact that confinement is not the only nonperturbative phe- 
nomenon that plays out in the QCD vacuum. The spontaneous breaking of 
approximate chiral symmetry is the other important player. Although we be- 
lieve that in principle the two phenomena are quite distinct, in the framework 
of QCD they are not easily separable. It is thus very important to have a 
simpler setting where one can understand the physics of confinement without 
having to worry about other complicated phenomena. 

Thankfully, we do know a model, or rather a class of models, which exhibit 
confinement but are free from other complications of QCD. Those are weakly 
interacting confining gauge theories in 2-1-1 dimensions. They were first studied 
by PolyakovB almost twenty years ago. By now we understand the confining 
physics of these theories very well and have several perspectives on it. These 
results have been confirmed by numerical methods of lattice gauge theoryfl 
Moreover in the last couple of years the finite temperature behavior of these 
theories has also been well understood, including the details of the deconfining 
phase transition. 

Our aim in this review is to discuss in detail, and from several comple- 
mentary points of view, the physics of these theories. We start in Sec. 2 with 
the description of the simplest of these models — the SU(2) gauge theory with 
the adjoint Higgs field — the Georgi-Glashow model. We describe Polyakov's 
derivation of confinement in this model based on the famous monopole plasma 
mechanism. In Sec. 3 we explain how confinement arises from the point of 
view of condensation of magnetic vortices. We show how the two pictures are 
closely related, but also point to some subtle differences between them and 
their complementary aspects. This has been discussed in detailEI in a previous 
volume in this series, and so our discussion will be brief. Section 4 is devoted 
to the discussion of Polyakov's confining string. We explain how the rigidity 
of the mathematical string arises from the absence of "stringy" degrees of free- 
dom in the ultraviolet. We also explain why the W bosons play the role of DO 
branes in the theory of the confining string. In Sec. 5 we show how the physics 
of confinement can be understood using Hamiltonian methods. We present a 
gauge invariant Hamiltonian variational calculation of the ground state wave 
function as well as of the potential between external sources. Section 6 is de- 
voted to a detailed discussion of the deconfinement phase transition and the 
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calculation of critical exponents. In Sec. 7 we discuss the properties of the 

vortex correlation functions and of the monopolc interactions in the high tem- 
perature phase. We show that the monopoles are linearly "confined" , while the 
vortex correlators decrease exponentially with distance. Section 8 discusses the 
same phase transition from the perspective of the confining string. It is shown 
that the would be "monopole binding" transition corresponds to the Hagedorn 
transition in the string language. The actual mechanism of the deconfining 
transition is the disintegration of the string due to proliferation of DO branes. 
Section 9 presents the generalization of various aspects of the preceding dis- 
cussion to the SU(A'') Higgs model. Finally, Sec. 10 contains some concluding 
remarks. 



2 The monopole plasma 

We start our discussion of confinement in 2+1 dimensions by recapitulating 
the original derivation of Polyakov. Consider the 2+1 dimensional SU(2) gauge 
theory with a scalar field in the adjoint representation, 



(1) 



We adopt the notation 



and 

D^h = df,h + g[A^,h]. 

The trace is normalized as 

trir'^T^) = 25"'' . 

We concentrate on the weakly coupled regime v ^ g"^. Since the coupling 
is weak one can try one's hand first at perturbation theory. Perturbatively 
the SU(2) gauge group is broken to U(l) by the large expectation value of the 
Higgs field. The photon associated with the unbroken subgroup is massless 
whereas the Higgs and the other two gauge bosons are heavy with the 
masses 

Mfj = 2Xv^ M^=g\,\ (2) 

Thus, perturbatively the theory behaves very much like electrodynamics with 
spin one charged matter. However, this is not the whole story. It is known 
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rigorousljfl that the weakly coupled regime in this model is not separated by a 
phase transition from the confining regime, which one naturally identifies with 
the vanishing Higgs expectation value. In the confining regime the spectrum of 
the theory is certainly massive, and the charged particles are linearly confined. 
Thus, there must be some nonperturbative effects that bring about these two 
changes even at weak coupling. The nonperturbative effects in question are 
the monopoles, as first discussed by PolyakovEl 

It is easy to see that the Euclidean equations of motion of the model 
have monopolc like solutions, 

D^Ff"" = g[h,D''h], Df,Dt'h= -X{h''h'' -v^)h. (3) 

The perturbative vacuum manifold is S0(3)/U(l) and this has a non-zero 
second homotopy group, n2(SO(3)/U(l)) = Z. Thus, the equations admit 
topologically nontrivial solutions with winding of the Higgs field at the space- 
time infinity. The "elementary" configuration of this type is the well-known 't 
Hooft-Polyakov monopolc, 

9A;ix) = \ea^.x^{l-c^{T), (4) 
where x"" = x"- /r. The non-Abelian field strength is 

9^-^ = - 1) + - e^.fcrx'')^', (5) 

where the prime denotes differentiation with respect to the 3D radial coordinate 
r. The boundary conditions satisfied by the various fields are h{0) = , 
h{oo) = V, 0(0) = 1 and </'(oo) = 0. The magnetic field far away from the 
origin is radial in the color space, and has a typical monopole-like behavior, 
dying away as In fact, the Abelian field strength 

F^, = h^'F;, ~ ^Eabch^iD^hfiD^hy (6) 

where 

N' 

is precisely that of a magnetic monopolc everywhere except a small region 
around the origin — the monopole core. The core size is of the order of the 
vector meson mass. At weak coupling the action of the monopole solution is 
very large 



2nA'L 



^monopolc- ^2 ^ j ' (7) 



(Mh\ 
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where 



Mw) 

is a slowly varying function such that 1 < e < 1.787 (further details can be 
found in Ref. and e(oo) = 1. Due to this large value of the action the 
contribution of the monopoles to the partition function is very small. Thus, 
typical configurations contain only very low density of monopoles and this 
dilute "monopole plasma" can be treated reliably in the fugacity expansion. 

The monopoles interact via the three-dimensional Coulomb potential. The 
classical action for a configuration that contains far separated monopoles and 
anti-monopoles is 

27r qaQb , 2ttMw ( Mh \ 2 /qn 
^=T2>.^? ^ + — ^MT^ >.9a (8) 



(9) 



where Qa = ±1. The monopole fugacity is determined as 



£, = exp 



2ttMw ^ f Mh \ 
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Even though this fugacity is tiny, the monopoles have a decisive effect on 
the low-energy properties of the theory. Since the W bosons are very heavy, one 
can neglect their effects with good accuracy. In this approximation, the parti- 
tion function of the theory factorizes into the partition function of the massless 
photon and that of the monopole plasma. The two factors can be combined 
into the pjartition function of a "dual photon" field. Following Polyakov's 
Derivationtl we start from the statistical sum Zm for Coulomb plasma with 

positive and iV_ negative charges with equal fugacities ^, which is given 

by 

^pias„.a= E [Ud\^aUd'r>'^''^''^^''-^l^^ (10) 

7V+,JV_-^ a=l fc=l +■ 

where fa and fb are coordinates for positive and negative charges respectively. 
The dimensionful constant fi is necessary to convert the summation over the 
monopole coordinates into the integration. It is nonuniversal, but clearly is 
related to the shortest distance scale at which the present approximation is 
valid, namely the size of the monopole core. In Polyakov's calculation fi is 
determined to be 

7/2 



/x^oc^. (11) 
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Using the fact that the CoTiIomb potential is the inverse of the Laplacian 
operator, this can be written asO 



^plasma = J i: y n ^'r^ n '^'^^^ nIinj. ^^^^ 



N+,N- " a=l b=l ' +■ 



exp(- /" dx^^-^{df,(j)f -Yi(l3{x)[{d{x - Xa) - S{x - Xb))] 



327T 

Here Zq is the partition function of a single massless field. In the dilute gas 
approximation, that is allowing at most one monopole or antimonopole at each 
spatial point, the summation over the coordinates is easily performed and one 
obtains for the partition function of the theory 



Z = J D(t)exp 
with the effective Lagrangian 



d'^xCc 



(13) 



^off = ^(^a'</')'-Cm'cos0. (14) 



Note that the factor Zq^ in Eq. (13) cancelled against the contribution of a 
free massless photon of the Georgi-Glashow model. 

One striking effect of the monopole plasma is explicit in this representation. 
The field 4> which in Eq. ( |l^ represents the photon, is not massless. Rather 
it has an exponentially small mass, determined by the monopole fugacity, 

(15) 

The other effect is confinement, or area law for the fundamental Wilson loop. 
Let us calculate the expectation value of the Abelian Wilson loop 



W{C) = (cxp (^ig J d^xB{x) 



(16) 



where i? is a component of the Abelian magnetic field (^) and S is a surface 
bounded by the curve C . 

To calculate the expectation value, we insert this exponential factor into 
the path integral for the partition function of the model and use the fact that 
for a monopole plasma configuration 



exp (^g j d^xB{x)^ — exp J d"^ 



xr]{x)p{x) . (17) 
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Vi-) = I TT^d'a.iy) (18) 



Here 

Is \x-y\' 

is the flux (normalized at 47r) through the surface E from a monopole located 
at the point and 



Pi^) ^^[6(X - Xa) - S{X - Xb)] 
a,b 

is the density of the magnetic charge in the plasma. As far as the summa- 
tion over the monopole coordinates is concerned this is exactly the same type 
of term as the linear term in in the exponential in Eq. (p^. Thus, the 
summation can again be performed, and one obtains 



WiC) = i / D<f>exp 



(fxCcsiv) 



(19) 



CMv)^7^idf.^?-^cosi^ + v)- (20) 



where 

For large contours C the path integral over can be calculated now in 
the steepest descent approximation. It is convenient to change variables to 
$ = (f>+r]. The integral then is determined by the solution of classical equations 
of motion 

92($-77) =M2sin($). (21) 

The nature of this solution is easy to understand. Let us for definiteness 
imagine that the surface E lies in the 1-2 plane. The Laplacian of 77 vanishes 
everywhere in space, except on the surface E, where it is equal to the transverse 
derivative of the delta function with the coefficient 2tt, 

d^f]{x) =2TTdJ^{x-J:). (22) 

Thus, the field $ satisfies the free sine-Gordon equation everywhere in space, 
except on the surface E. The presence of rj imposes the boundary condition 
that this particular solution must be discontinuous across E with a disconti- 
nuity of 27r. Such a solution is easily constructed from a domain wall solution 
of the sine-Gordon theory. At negative z the solution $ is just the domain 
wall, so that at z = 0-1- it reaches the value tt. It then jumps discontinuously 
to — TT on the other side of the surface, $(2; = 0—) ~ — vr and rises back to 
zero at large positive z again following the domain wall solution. The action 
of such a solution is obviously just equal to the action of a domain wall in the 
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sine-Gordon model. It is proportional to the total area of the surface S, and 
the coefficient parametrically is given by g^Af^. Thus, we find 

W(C) - exp {-(tS) (23) 

where S is the minimal area and cr ~ g^M^. The Wilson loop therefore has an 
area law and the theory is confining. 

Thus, the monopole plasma accounts nicely for the finite mass of the pho- 
ton as well as for the area law of the Wilson loop. 



3 The vortex condensation 

A complementary view of confinement in this theory is the spontaneous breakn-, 
ing of the magnetic Z2 symmetry due to the condensation of magnetic vortices™ 
This approach has been reviewed in Ref. ^ and therefore we will be brief here. 

The Georgi-Glashow model has a global Z2 symmetry. This symmetry is 
not recognizable as a Noether type symmetry in the Lagrangian (|l|) , but rather 
is akin to "topological" symmetry. It is generated by the total magnetic flux. 
Ignoring for the moment monopole effects, the Abelian magnetic field (^) is a 
conserved current 

= \^^..xF,x. (24) 

However, due to the monopole contributions (see Ref. p^ , this symmetry 
is anomalous, and only the Z2 subgroup is. conserved. The non-anomalous 
Z2 magnetic symmetry group is generatedEJ by the large spatial Wilson loop 
which encloses the system, 

W{C 00) = exp (^-^ j (fxB{x)^ . (25) 

The order parameter for this Z2 transformation is the operator that createsB0 
an elementary magnetic vortex of fiux 2T:/g, 

V{x)^exp / e,jh''E''{x). (26) 

9 Jc 

This operator is local, gauge invariant and a Lorentz scalar. Together with the 
spatial Wilson loop it forms the order-disorder algebra, 

W{C ^ oo)Vix) ^ ~V{x)W{C ^ 00) . (27) 

The vortexj-jOperator condenses in the vacuum and its expectation value is 
determinedcl by the gauge coupling constant (V)'^ = g^/Sir'^. 
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Since the anomalous effects tliat break the U(l) magnetic symmetry down 
to Zi at weak coupling are small, one can construct the low-energy effective 
theory Lagrangian in terms of the relevant order parameter. It is given by the 
Z2 invariant theory of the vortex field V ^ 

8tt^ 4 



c,s = d^vd^'v* -\[vv* -^j +—{v' + (yy}. (28) 

E 



The vortex self-coupling in this effective Lagrangian can be determined as 

A = 29 

9 



Note that as a low-energy Lagrangian, Eq. 1)^) is indeed consistent with 
Eq. (p^. At weak gauge coupling, the quartic coupling A is very large. In 
this nonlinear cr-model limit the radius of the field V is therefore frozen to its 
expectation value. The only relevant degree of freedom is the phase 

Vix) = -^e^^. (30) 

Substituting this into Eq. ( |2^ ) one indeed obtains precisely Eq. (|l^) with the 
identification (p = 2x- Thus, Polyakov's dual photon is nothing but the pseudo- 
Goldstone boson of the Lagrangian (pF 



This relation can be inferred directly from the definition of V in Eq. (26) 
and the monopole plasma partition function (|l3|). Equation (13) relates the 
field (j) with the monopole density as 

£,m-^P-id,P^- (31) 

The factor i has disappeared in the second equality since the field strength 
is now in Minkowski rather than Euclidean space. The factor ^ is due to 
the fact that the minimal magnetic monopole charge is twice the fundamental 
value allowed by the Dirac quantization condition. This immediately gives 

= — = — / dz^F^^iz). (32) 



3 92 g 



Xq 



The integral on the right-hand side does not depend (to leading order in fu- 
gacity) on the choice of the integration curve. Choosing it to be at equal time 
and comparing to Eq. (|2^), we get 

= 2x. (33) 
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Thus, the non- vanishing photon mass is a natural consequence of the fact that 
the magnetic symmetry is discrete and not continuous, and thus the photon is 
a pseudo-Goldstone boson. 

The calculation of the spatial Wilson loop in this framework closely follows 
the calculation in the previous subsection. For details we refer the reader to 
Ref. |. 

Thus, as far as the description of the low-energy phenomena is concerned, 
the vortex approach is equivalent to the monopole plasma approach. The one 
conspicuous difference between them is that in the vortex low-energy approach 
the field x is treated as a phase, while in the monopole plasma context (p is an 
ordinary real field which takes values on a straight line. As a consequence the 
Lagrangian ( |28| ) allows configurations with quantized discontinuities of the field 
X- In particular, the vortex configurations where the field x winds around some 
point are allowed in the description of Eq. (psj). In such a vortex configuration 
the field x is discontinuous along an infinite cut. In the framework of the 
effective theory, Eq. (p^), such a quantized discontinuity does not contribute 
to the energy. Polyakov's effective Lagrangian (p^, on the other hand, does 
not allow such vortex configurations. Since the field (j) in it is not treated as 
a phase, the cut does contribute to the kinetic energy term in Eq. (p^, and 
this contribution is both ultraviolet and infrared divergent. The ultraviolet 
divergence is of no importance by itself, since the theory is defined with the 
intrinsic ultraviolet cutoff Mw- However, the infrared divergence indisputably 
puts the description of such vortex- like configurations beyond Eq. (p^. 

As it turns out this difference can sometimes be quite significant. The 
point is that these vortex-like configurations of the field V are the low-energy 
incarnation of the charged states of the theory, the W bosons. 

As explained in Ref. ^, the electric current is identified with the topological 
current in Eq. (|2^), 

'':Jf. = ief..xd,{V*d)y). (34) 



TT 



Thus, charged states carry unit winding number of %. The W'^ [W") boson 
is a state with positive (negative) unit winding of x- 

Thus, the physics of the charged states cannot be discussed within Eq. (Q), 
but can be considered within Eq. (|2|). Of course since the charged states 
are confined in this theory, the energy of an isolated charged state is linearly 
infrared divergent. The mechanism of this divergence, as well as the formation 
of the electric fiux tube, is discussed in detail in Ref. ^. As we will see in later 
sections, proper treatment of the charged states is crucial for the understanding 
of the dynamics of the deconfinement phase transition and thus it is worthwhile 
to be careful on this point. We note that this is not the only instance in which 
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the charged states play an important role. For example, as discussed in Rcf. 
O, the presence of these states leads to the breaking of the string of charge 



two, which in Polyakov's effective Lagrangian (14) is strictly stable. 

In fact, to properly account for the large mass of the charged particles, the 
Lagrangian ( p8[ ) must be slightly modified, since as it stands it underestimates 
their energy. The core energy of the soliton of V should be equal to the mass 
of the charged vector boson Mw, since this is the lightest charged excitation 
in the theory. On the other hand, the core mass of a soliton in the Lagrangian 
( p8| ) is given by the ultraviolet contribution of the Coulomb potential. With 
the cutoff of order Mw, this is of order In Mw/g'^- This is not surprising, 
since the mass of the W boson indeed comes from the distances of the order 
of its Compton wavelength, and is thus well inside the ultraviolet cutoff of 
the effective theory. The situation can be improved by adding to the effective 
Lagrangian a higher derivative term of the Skyrme type 

dL = A{ef,,xd,V*dxV)^, (35) 

with 

A(xl-^. (36) 

The extra term Eq. (|3^) does not affect the photon mass nor the value of the 
string tension. The density of this action does not vanish only at the points 
where the phase of V is singular — that is the points of winding. For a closed 
curve C of length L which carries the winding, the contribution of this extra 
term to the action is MwL, which is precisely the action of the world line of a 
massive particle of mass Mw ■ 

Thus, Polyakov's monopole plasma picture and 't Hooft's magnetic vortex 
condensate picture fit very neatly together and provide a fairly exhaustive 
understanding of confinement in the Georgi-Glashow model. Of course, both 
monopoles and vortices are not physical excitations of the Georgi-Glashow 
theory. What one would like to see is the confining string, the tube of electric 
flux. It is believed that in 3+1 dimensional QCD at large Nc these strings 
interact weakly and perhaps are described by a free string theory. Although 
one does not expect the same to be true in 2+1 dimensions at Nc — 2, the 
handle we have on the present model gives us a unique possibility to try and 
derive the relevant string theory. In the next section, therefore, we explain 
how one can derive in the long wavelength approximation the theory of the 
confining string. This theory turns out to have some quite peculiar features, 
which we will discuss. 
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4 Confining strings 

In 1996 Polyakov proposedQ the derivation of the so-caUed confining string 
action expHcitly fi'om compact QED, which is equivalent to the Georgi-Glashow 
model in the low-energy approximation of Eq. (p^). This aciim was studied 
further in Ref. |l^ and it was found to describe a rigid stringE3llZl 

We will first discuss, the original derivation due to Polyakov which uses 
the Kalb-Ramond fieldliS We will also give a more intuitive and simple deriva- 
tion which uses the direct correspondence between field configurations and the 
summation over closed string world sheets. An advantage of this derivation is 
that it makes explicit an important and quite unusual property of the confining 
string, namely that the fluctuations of this string with large momenta (larger 
than the inverse thickness of the string) do not cost energy. This is a dynam- 
ical extension of the so called zigzag symmetry introduced by PolyakovEJ We 
also show that the heavy charged particles of the Georgi-Glashow model {W^ 
bosons) appear as a certain kind of branes in the string description. 

Our starting point will be an expression for the Wilson loop discussed 
earlier, 



WiC) = i y" i?0exp ~J d'xC^siv) 



(37) 



where 

CM = -i^id,<pf - cos(</) + (38) 

Note that the effective theory at 77 = has more than one vacuum state. In 
particular, = 27rn for any integer n is the classical ground state. Therefore, 
the classical equations of motion have wall-like solutions — where the two 
regions of space, say with cf) = Q and (f> — 2Tr are separated by a domain wall. 
The action density per unit area of such a domain wall is parametrically 

a cx g^M. (39) 

As we discussed in Sec. 2, the fundamental Wilson loop, when inserted into 
the path integral, induces such a domain wall solution, with the result 

{W{C)) = exp{-aS}, (40) 

where 5* is the minimal area subtended by the contour C. Thus, the string 
tension is precisely equal to the wall tension of the domain wall. This domain 
wall is identified with the world sheet of the confining QCD string. 
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The Kalb-Ramond way 

In order to find a string representation for W{C) one can rewrite W{C) in a 
different form, introducing a new auxiliary Kalb-Ramond field B^^, 



WiC)= J DB^,Dq>exp[-S{q^,B)] 



(41) 



with the action 



S{c^,B) = J 



16^2 M^g\ 

■Bt^ + ^^t^iypOi,(t>Byp + - COS( 



r 



167r2 



+ i 



I B/^i^da^^. 



(42) 



Integration over B brings one back to the initial expression for Wilson loop 
Eq. (|3^). We may choose, however, to integrate first over (p. In that case we 
obtain an effective action for the massive B field, 



W{C)=JdBcxp j dx(^-^Bl^ + J{H)^ exp*^ 



Bda, 



(43) 



where H = dB is the field strength and the function H is given by (for details 
see Ref . |l|) 



H arcsin 



H 



(44) 



The multivaluedness of this action reflects the periodicity in (j) mentioned ear- 
lier. Expanding in powers of B and limiting ourselves to quadratic terms, the 
effective action takes the form. 



d^x 



Bi 



const _ 
~M2 



S{B)- 
For large loops one gets 

I^(C)~exp(-F(C, S))) 
with the following local expansion for the string action F: 

F = cig^M J d^^^ + C2gH4-' J d^a^t^,f^ + 



where 



gab = daXdbX, t^u = e° daX^dbXi, 



(45) 

(46) 

(47) 
(48) 
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This is the action of the rigid string which also can be written as (see Ref. 
and references therein) 

F = J d'^^Vgg°'^DaX (fj - sD^ + ..) DbX (49) 

where g and Da respectively are the determinant and the covariant derivative 
with respect to the induced metric gab = daxdtx on the embedded world sheet 
a;(^i,^2)- Note that the rigidity is controlled by the second term and one can 
show that the stiffness s is parametrically of order 

so.^. (50) 

The stiffness s in this expression is actually negative, but the system is 
stable due to higher order terms. Note that since a oc g'^M ^ M^, the rigidity 
is very large, and is effective on the distance scales much larger parametrically 
than the "natural" string scale given by the string tension. 

The action (^^ explicitly depends on the surface E , while originally it was 
introduced as an unphysical object. The reason for this is the multivaluedness 
of the full action S{B). As was shown in Ref. |lj the surface independence of 
{W{C)) is restored after summation over all possible surfaces. 



4-2 The confining string for pedestrians 

Let us now derive the action of a fluctuating string in a More direct wayEl We 
work with the Polyakov effective Lagrangian (|lj) . Since the string world sheet 
is identified with the domain wall in the effective action, we will integrate in the 
partition function over all degrees of freedom apart from those that mark the 
position of the domain wall. To do this let us split the field (f> into a continuous 
and discrete part 

(l){x) = 4>{x) + 2TT'q{x) (51) 
where the field (j) is continuous but bounded within the vicinity of one "vacuum" 

- TT < (^(x) < TT (52) 

and the field rj is integer valued. Clearly, whenever rj{x) ^ 0, the field (j) is 
necessarily not in the vicinity of the trivial vacuum (j) = Q. Thus, the points 
in space where rj does not vanish mark in a very real sense the location of a 
domain wall between two adjacent vacua. The partition function Z is given 
by the path integral 

i9^((^ + 27177)9^(0 + 2^7?) 
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+ M^(l-cos(^) 



(53) 



At weak coupling the field <^ can be integrated out in the steepest descent 
approximation. This means that the equations of motion for cj) have to be 
solved in the presence of the external source 77, keeping in mind that (j) only 
takes values between — tt and tt. Let us assume for this solution that the 
surfaces along which r} does not vanish are few and far between. Also, we 
will limit the possible values of 77 to 0, 1, —1 corresponding to the "dilute gas" 
approximation. In fact, as we will see below, this exhausts all physically allowed 
values of t}. Then the qualitative structure of the solution for is clear. In the 
bulk, where rj vanishes, the field (f) satisfies normal classical equations. When 
crossing a surface S on which 77 = 1, the field (p must jump by 2-k in order to 
cancel the contribution of r] to the kinetic term, since otherwise the action is 
ultraviolet divergent. Thus, the solution is that of a "broken wall". Far from 
the surface the field ^ approaches its vacuum value (j) = Q. Approaching S it 
increases to tt, then it jumps to — tt on the other side of the surface, and again 
approaches zero far from the surface. The profile of the solution for smooth S 
(in this case a plane) is depicted in Fig. 1. Since outside the surface the field 
(f) solves classical equations of motion, clearly the profile of (j) is precisely the 
profile of the domain wall we discussed above. As noted above, the purpose 
of the discontinuity across the surface S is to cancel the would be ultraviolet 
divergent contribution of r} to the kinetic term in the action. Thus, the action 
of our solution is precisely 

S[„„s = ./.v|^ai, (54) 

where t are the coordinates on the surface S, and are the coordinates in 
the three-dimensional space-time. This is precisely the Nambu-Gotto action 
of a free string. 

Of course the confining strings in the Georgi-Glashow model are not free. 
The thickness of the region in Fig. 1 in which the field (j) is different from 
zero, is clearly of order of the inverse photon mass M~^. Thus, when two 
surfaces come within distances of this order they start to interact. In principle 
this interaction can be calculated by just finding the classical solution for (j) in 
this situation. Now however we want to discuss one particular property of the 
confining strings — their rigidity. 

4-3 The string is soft ... and therefore rigid! 

The Nambu-Gotto action we have derived in the previous subsection is of 
course only the long wavelength approximation to the action of the confining 
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Figure 1: Schematic representation of the solution profile of the ^ field as a function of the 
s-coordinate. 



17 



Figure 2: The solid line represents the actual thin string. The dashed lines denote the 
contours of the effective thick string. 

string. It is only valid for string world sheets which are smooth on the scale 
of the inverse photon mass. Expansion in powers of the derivatives can in 
principle be performed and it will give corrections to the action on scales 
comparable to M^^. However, for our confining string, strange things happen 
in the ultraviolet. Physically, the situation is quite peculiar, since the action of 
the string has absolutely no sensitivity to changes of the world sheet on short 
distance scales. This should be obvious from our derivation of the string action. 
Suppose that, rather than taking 77 = 1 on an absolutely smooth surface, we 
make the surface look the same on the scale Af~^ but add to it some wiggles 
on a much shorter distance scale d, as in Fig. 2. To calculate the action we 
now have to solve the classical equation for the field (f> with the new boundary 
condition — the surface of the discontinuity is wiggly. This new boundary 
condition changes the profile of the classical solution only within the thickness 
d of the old surface. However, since the action of the classical solution is in no 
way dominated by the region of space close to the discontinuity, the action of 
the new solution will be the same as that of the old solution with the accuracy 
dM^. Thus, all string configurations which differ from each other only on small 
resolution scales, d ^ Af.^^, have to this accuracy the same energy! The string 
is therefore extremely soft, in the sense that it tolerates any number of wiggles 
on short distance scale without cost in energy. Fig. 3. In particular, since the 
string tension for our string is much greater than the square of the photon 
mass, (j/M?^ — /M^ 3> 1, fluctuations on the scale of the string tension are 
not penalized at all. 

This independence of the action on short wavelength fluctuation is a dy- 
namical manifestation of the so called zigzag symmetry introduced by PolyakovE 
Indeed, Polyakov notes that if a segment of a string goes back on itself, phys- 
ically nothing has happened and so such a zigzag should not cost any action. 
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This situation is an extreme example of a wiggle we have just discussed — a 

wiggle with infinite momentum. What happens physically is that not only do 
infinite momentum modes cost nothing, but so do finite but large momentum 
string modes. 

The confining string is therefore very different from a weakly interacting 
string usually considered in the string theory. In the weakly interacting string, 
momentum modes with momentum of order of the square root of the string 
tension carry energy which is of the same order. In the confining string on 
the other hand, these momentum modes do not carry energy at all. Thus, 
we do not expect the spectrum of the confining string to contain states with 
large spatial momentum (heavy states with low angular momentum). In fact, 
the states with high angular momentum may very well be absent too, due to 
the fact that a long string can decay into a — W~ pair. Therefore, it is 
likely that most of the string states will not appear in the spectrum of the 
Georgi-Glashow model. Thus, our confining string theory is quite peculiar: it 
is a string theory without stringy states in the spectrum. 

In a somewhat perverse way, this softness of the "mathematical" string 
leads to rigidity of the physical string. The point is the following. As is obvious 
from the previous discussion, the high momentum modes are in a sense "gauge" 
degrees of freedom. The action of any string configuration does not depend on 
them. Consider a calculation of some physical quantity O in the string path 
integral. If O itself does not depend on the string high momentum modes, the 
integration over these modes factors out. If, on the other hand, O does depend 
on them, then its average vanishes since their fiuctuation in the path integral 
is completely random. This is absolutely analogous to the situation in gauge 
theories, where all observables must be gauge invariant, and for calculation of 
those the integration over the gauge modes always factors out, independently 
of the observable under consideration. 

Thus, for all practical purposes, we should just exclude the high momen- 
tum string modes from consideration altogether. For example, when calcu- 
lating the entropy of our string, we should only take into account the states 
which are different on the coarse grained level with the coarse graining scale of 
order M~^. This means that our string is intrinsically "thick". If we still want 
to describe this situation in terms of a mathematical string with the string 
tension a, we must make the string rigid so that any bend on the scale smaller 
than is suppressed. Thus, such a string theory must have a curvature 

term with the coefficient of order cr/M^ 



(55 oc (curvature). (55) 
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This is indeed what happens in the truncated confining string action (p9|). 

4-4 The bosons — the branes of the confining string 

Interestingly enough, the charged states of the Georgi-Glashow model have a 
natural place in the confining string framework. To represent a W boson in 



the string language, let us consider a path integral Eq. ( |14[ ) in the presence 
of one soliton. As we discussed above, within the effective Lagrangian (14) 
proper, solitonic configurations do not exist. However, they can be forced onto 
the system with the help of the external current. To create a winding state 
with world line C we should insert a source in the path integral which forces a 
unit winding on the field x- In other words, the field (f) would have to change 
by An when going around C. The relevant partition function is 

Z[C] - J D(j)e^p {-MwL - j]) (56) 



where 



lid^cj) - 2nj^){d^(t> - + M2(1 - cos< 



(57) 



.2 

Here the "external current" is 

j^ix) = nl{x)S{x e Si) + nl{x)S{x € S2) (58) 

where 5*1 and S2 are two surfaces which both terminate on the curve C, and 
the unit vectors and are the normal vectors to these surfaces. The shape 
of the surfaces on which the current does not vanish is illustrated in Fig. 4. 
The insertion of forces the field (f> to jump across the surface 5*1 by 2n 
and again jump by 2tt across S2 in order to cancel the, otherwise ultraviolet 
divergent contribution of j^. Thus, when going around C the field (j) changes 
by Att and therefore C is the world line of the boson. The position of 
the surfaces and 52 is arbitrary as long as they both terminate on C. In 
particular, they could coincide, but here we prefer to use a more general form 
with non-coinciding surfaces. 

As before, splitting the field variable 4> into (j), and 77 we see that the 
presence of just shifts the variable 77 by unity on the two surfaces Si and 
^2. The integration over 4> at fixed rj is performed in exactly the same way as 
before. The only difference now is that for any given 77 one has two extra string 
world sheets along S'l and S'2. After the integration over 77, the result will not 
depend on the exact position of Si and S'2. However, at fixed 77 the two surfaces 
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introduced in Eq. ( pq ) specify the positions of the two extra world sheets in 
the confining string path integral. We thus see that the field-theoretical path 
integral in the presence of a W'^ boson, in the string representation, is given 
by the sum over surfaces in the presence of a — brane, which serves as a 
source for two extra string world sheets. 

Let us stress here the similarity and differences between W particles viewed 
as 0-branes, and DO branes in string theory. Usually DO branes are thought of 
as infinitely heavy. The situation in the Georgi-Glashow model is very similar 
in this respect. They are not infinitely heavy, but very heavy indeed since the 
mass of W is large on all scales relevant to zero temperature physics. However, 
this analogy with D branes in superstring theory is not exact. The mass of 
a DO brane is proportional to 1/gstr whereas the mass of the W boson (our 
brane) is proportional to g. The fact that it is very heavy in the weak 
coupling limit is due to the large value of the Higgs condensate. It may be 
more natural to think of the as of a stretched F-string, which only after a 
duality transformation becomes a D-object on which the confining string can 
end. For the sake of brevity, in the following we will disregard this subtlety 
and will continue to refer to the W bosons as branes. 

The contribution of the brane to the partition function is suppressed by 
a very small factor exp{— Mvi/^}, and vanishes for an infinitely large system. 
As we will see in the later sections however, the situation changes dramatically 
at finite temperature, where one dimension of the system has finite extent. 

Incidentally, going back to our definition of the vortex field V , we see that 
in between the two world sheets the value of V is negative, while outside it 
is positive. Thus, we have created a domain of the second vacuum of V in 
between 5*1 and 5*2. It may be easier to visualize the situation with both 
and W" present. In this case the surfaces 5*1 and 5*2 terminate on one side 
on the world line of , and on the other side on the world line of W~ (see 
Fig. 5). They are thus boundaries of a closed domain of the second vacuum 
of the field V . In the infrared therefore our strings are nothing but the Ising 
domain walls, and the pair of DO branes creates an Ising domain. 

Note that in physical terms there are only two distinct vacua in the model: 
(y) = 1 and (y) = — 1. Thus, having two domain walls is the same as having 
a wall and an anti-wall, and if they coincide spatially such a configuration is 
equivalent to the vacuum. A configuration of rj with 77 = 2 on a closed surface 
is physically equivalent to vacuum. Therefore, the values that rj is allowed to 
take are limited to 0, 1 and —1. 

We close this section by noting, that the reason we have two string world 
sheets terminating on a DO brane is that the field theory in question only 
has fields in the adjoint representation of SU(2). We can imagine adding 
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heavy fundamental particles to the model. We would then also have allowed 
configurations of one string world sheet terminating on a DO brane. These 
DO branes would however be physically different and would have a different 
mass and therefore a different weight in the path integral. In field-theoretic 
terms, presence of the fundamental charges drastically changes the properties 
of the Z2 joiagnetic symmetry, turning it into a local, rather than a global, 
symmetryEj 

5 The Hamiltonian approach to compact U(l) 

Our discussion so far has been in the context of the Lagrangian path inte- 
gral formalism. Often, Hamiltonian approaches tend to be useful and more 
instructive. In this section we will therefore show how the same confinement 
phenomenon manifests itself in the canonical Hamiltonian analysis. 

To make life a little easier in this section we will not study the SU(2) 
Georgi-Glashow model, but rather compact U(l) theory. This theory is also 
confining, and at low energies is indistinguishable from the Georgi-Glashow 
model. Since we will not at this point ask questions about the charged sector 
of the theory, compact U(l) is perfectly adequate for our purposes. In this 
section we follow Refs. ^ and 

We start therefore by setting up the Hamiltonian description of U(l) com- 
pact QED. The first question to be settled, is what is the Hilbert space of 
admissible states. Firstly, it is clear that the Gauss' law should be imple- 
mented, and thus all the physical states should satisfy 



There is a crucial difference between the Gauss' law in the compact theory 
and in the noncompact one. In the noncompact theory equation (^9|) should 
be satisfied only for regular functions A. For example, the operator 



which has the form of (|9|) with the function A proportional to the planar 
angle 0, i.e. A = -Q{x), does not act trivially on physical states. In fact, this 
operator creates point-like magnetic vortices with magnetic flux lix j g just like 
Eq. (p6|). It therefore changes the physical state on which it is acting. 

In the compact theory the situation in this respect is quite different. Point 
like vortices with quantized magnetic flux 2iTn/g cannot be detected by any 
measurement. In the Euclidean path integral formalism of Ref. this is the 




(59) 




(60) 
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statement that the Dirac string of the monopole is unobservable and does not 
cost any (Euchdean) energy. In the Hamiltonian formahsm this translates 
into the requirement that the creation operator of a point-hke vortex must 
be indistinguishable from the unit operator. In other words, the operator 
( |60| ) generates a transformation which belongs to the compact gauge group, 
and should therefore act trivially on all physical states. Equation ( |5^ should 
therefore be satisfied also for these operators. 

Accordingly, the Hamiltonian of the compact theory also must be invariant 
under these transformations. The magnetic field, defined as -B = CijdiAj, on 
the other hand does not commute with V{x) (cf. Ref. [2^ ), 

V\x)B{y)V{x) = B{y) + jS^x - y). (61) 

The Hamiltonian should therefore contain not but rather a periodic function 
of B. We will therefore choose our Hamiltonian to be 

^=^«'E^--^Eco«5fl'i?n. (62) 

Since we will need in the following an explicit ultraviolet regulator, we use lat- 
tice notations. Here a is the lattice spacing, and the sums are respectively over 
the links and plaquettes of the two-dimensional spatial lattice. The coefficients 
of the two terms in the Hamiltonian are chosen so that in the weak coupling 
limit, upon formal expansion to lowest order in , the Hamiltonian reduces to 
the standard free.Hamiltonian of 2-1-1 dimensional electrodynamics. Following 
Polyakov's bookP we work in the weakly coupled regime. Since the coupling 
constant in 2-1-1 dimensions has dimension of mass, weak coupling means 
that the following dimensionless ratio is small 

5^0 <1. (63) 

Our aim now is to find the vacuum wave functional of this theory. Exact 
solution of this problem is not feasible. However, at weak coupling we can use 
the variational approximation. 



5.1 The variational ansatz 



For a weakly coupled theory one expects the vacuum wave functional (VWF) 
to be not too different from the vacuum of a free theory. Since the VWF of 
free (noncompact) electrodynamics is Gaussian in the field basis. 



'(/'o[A] = exp 



i^^r.G-i(r-s)^ 

r,s 



(64) 
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the Gaussian variational approach in this case should give a good approxima- 
tion. An important caveat, however is that the ground state WF should be 
gauge invariant under the full compact gauge group. As a result it turns out 
that one can not take just a Gaussian in A, since this will not preserve gauge 
invariance. The simplest generalization of the Gaussian ansatz which we use 
along the lines of Ref. is to project a Gaussian WF into the gauge invari- 
ant subspace of the Hilbert space. The projection has to be performed with 
respect to the full compact gauge group of Eq. (59). 

To facilitate this, define a vortex field A^^ that satisfies (we suppress the 
lattice spacing a henceforth) 



"5 



V • A^ = 



(65) 



where n' is a plaquette on the lattice, or, equivalently, a site of the dual 
lattice. This is the vector potential corresponding to a magnetic field that is 
zero everywhere except at n' = 0, where it takes the value The explicit 
solution of ( |65| ) is 

AYin) = -^e., ( ^6^,,o] (66) 
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The compact gauge invariance requires that the variational wave function ip[A] 
be invariant under shifts A — > A -f A^. This of course is consistent with the 
periodicity of H under B —t B + We also demand noncompact gauge in- 
variance of the wave function. Hence we define a field, shifted by a noncompact 
gauge transformation 0„ and by a vortex distribution rrij^i , 



A(f'") = A„ - (V(^)n - ^m„,A^'(n - n') 



(67) 



or, for short, 

^(0,m) _ (68) 

We choose the gauge invariant and periodic trial wave function as 



^[A]= /[d0n]expU^4f")G-nr-s)4r^ 



Under a gauge transformation. 



(69) 



V'lA-f VA] = ^[A] 



(70) 
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since A can be absorbed in a shift in (j). The simple rotational structure of 
Gij = SijG that appears in the variational wave function ( |69| ) is consistent 
with perturbation theory. We also take G{x) to be a real function. 

Our task now is to calculate the expectation value of the Hamiltonian in 
this state, and to minimize it with respect to the variational function G. We 
start by considering the normalization of the wave function. 



5.2 Normalization integral 
The normalization of \ip) is 



X exp 



1 



{m,m'} ' 



exp 



}_^{4>' ,m')Q~l j^{4>' ,m') 



(71) 



We shift A by V(/)' + A • m', and absorb the shift into and m, giving 



z = j2 jmwe- 

{m} 



e 2 



i,AG~'A 



(72) 



Now we combine the exponents according to 

j^{<t>;m) Q-l j^{(j>,'m) yi^Q-l _ 2^('^'/2,™/2) (^-1^(0/2, m/2) 

1 



+ -S{(l>,m)G-^S{(t>,m). 



where 



S = V0 + A*^ - m . 



(73) 



(74) 



The first term in S has zero curl while the second is divergenceless. Fur- 
thermore, G~^ is translation invariant and proportional to the unit matrix. 
Thus, SG^^S has no cross terms between m and (f>. We now shift A by 
V0/2 + A^ • m/2, and all the fields decouple. We have then 



where 



Z Zj\Z^Zy 

Za = det ttG 



(75) 
(76) 



[d(l)]e- 
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= (dct47r;^G 



1/2 



r2 X! ™r'£'(r' - s')ms/ 



(77) 
(78) 



Z„ = ^ cxp 
{'"„'} 

Here is the vortex partition function, with the vortex-vortex interaction D 
given by 



D(r' - s') = g^Y^ A^(r - r') • G-^[v - s)A^(s - s') , 



or 



D 



We can spHt off the r' = s' terms in (^Sh and write 



H exp 

{™n' } 



n 



where we have defined the vortex fugacity 

2 — e 49-! 



(79) 
(80) 

(81) 

(82) 



In the interest of clarity, we adopt henceforth a continuum notation, indi- 
cating the uhraviolet cutoff a only where necessary. Moreover, since the inter- 
esting physics comes from the infrared, lattice effects can be approximated by 
a momentum-space cutoff A = a~^, which can simplify formulas further. 

The variational function G appears both in the vortex-vortex potential 
and in the vortex fugacity. We expect the ultraviolet behavior of G to be the 
same as in the free theory, viz., 



so 



D{Q) 



G-\k) - k , 



(Pk 47r2 



{2n 



,2 1,2 



G'\k) 



27rA 



(83) 

(84) 
(85) 



and thus 



z - e ^ (86) 
In the weak coupling region we have z <C 1, justifying a restriction to m = 0, ±1 



in (|78|) and (|8^ 
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5.3 Expectation values 

We calculate correlation functions of m via a duality transformationi We add 
an i J • m term to the exponent in (fsl) and use the formula 



e 4g2 ™ ™ _ const I [dx\ 



(87) 



to obtain 



Noting that0 



Zv= j [dx] e-s"^-'"'"-^ + 2cos(x„ + J„) 



cos(x + J) = (cosx)o :cos(x + J):= z :cos(x + J): 



(89) 



we have 

Z,, 



[][l + 2z:cos(x + J):] 



Dxexp -g^xD ^X + 2z j d^x : cos {x{x) ^ J [x)^ : 

in continuum notation. Correspondingly^ 

{m{x)m{v)) = 2g^D-\x ~ y) ~ 4/ {D-\(x) D-\{v)) . 



(90) 



(91) 



The propagator of x is easily calculated. To first order in z, the only contri- 
bution comes from the tadpole diagrams, which have already been subtracted 
in Therefore, 

d^xe^-{x{x)xm 



D{k) D^ik) 



The correlation function of the vortex density is then 



(92) 



K{k) = J d'^x {m{x)m{Q)) = 2z + 0{z^) , (93) 
which in this approximation does not depend on momentum. 

"T hp n nrmal ordering is performed relative to the free theory defined by the quadratic action 
in H). 
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Now we are ready to calculate the expectation value of the Hamiltonian 
(|6^). Using the definition ( |69| ) we obtain 



1 E^d^x 


)- 




(2^)2 


\k) 


9' J 


(2^)2 




2^2 

~^ 

9^ 



d^k 

d^k 
(2^ 



k-^G-^{k)K{k) 



(94) 
(95) 



The magnetic part is easily calculated since it has an exponential form and, 
therefore, with our trial wave function leads to a simple Gaussian integral. We 
find 

d^k 



(e™9Bn^ =exp 



1 2 2 

-—nq 
4 ^ 



{2n) 



;k^G{k) 



(96) 



The second factor is due to the vortices, and is different from unity only for 
odd values of n. Using ( pO| ) we find easily that ^e*'^™) — e^*^. Expanding to 
leading order in and z, we get 



9' 



■ cos gB 



d'^k ,2^,,^ 4 
k^G{k) + —z , 



(27r) 
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(97) 



where we have dropped an additive constant. Finally, the expression for the 
variational vacuum expectation value of the energy is 



d^k 
(27r)2 



G-^{k) + k^G{k) 



47r" 



z k-^G-^ik) 



(98) 



5.4 Determination of the ground state 

The expression ( p8|) has to be minimized with respect to G. From equation 
(H) and (11) we find 

^ = 4^^"^^"^(^)^ 
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The variation of (^8|) gives|j 

9 J (27r 



Equation (100) has the solution 



where 



47r4 



d'k 



fc2 



(100) 



(101) 



(102) 



The main contribution to the integral in the gap equation (102) comes from 
momenta k^ 3> nn? ■ For these momenta k'^G~^{k) = 1. Therefore, we see that 
(102) has a nontrivial sohition. Using Eqs. (p3), (p4), and (101) we obtain 



An" 



■exp - — 



J (27r)2^^ 



_fk_ 
(2^ 



fc2 



which for g^a = /A ^ 1 can be simplified to [cf. (|8^)] 



^ 2(^'-4)A4 
Att — exp 

9^ 



ttA 



(103) 



(104) 



where we have restored the ultraviolet cutoff dependence exphcitly. The re- 
sulting m is the mass gap of the theory, in the sense that it is the inverse 
of the spatial correlation length. Calculating, for example, the propagator of 
magnetic field, we find 



and at large distances (neglecting power-like prefactors), 

y^G{x) = - / + m2)i/2e'k.x _ g-,n. _ 

J V'T^Y 

This dynamically generated mass is Polyakov's resulti Thus, we recover in 
the Hamiltonian approach the first important result — a finite mass gap m. 



(105) 



(106) 



We have dropped a term 



-zk G (fc) from the right-hand side of (100) is smaller by 



a factor of than the term retained, assuming G ~ fc ^ at large k. 
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Let us note that a mass gap m in a Hamiltonian formalism (104) has a dif- 
ferent preexponential factor in comparison with a mass gap M-y obtained in 
a path integral approach (^5|). The reason is very simple. Contrary to uni- 
versal exponential factor the preexponential factor depends on details of UV 
regularizations which are different in Lagrangian and hamiltonian formalisms. 

5.5 Spatial Wilson loops 

We also want to see whether the charges are confined in our best variational 
state. The simplest quantity that is related to confinement is the expectation 
value of the Wilson loop. Therefore, we will calculate it in our ground state, 

Wi[C] = ^exp (^Ig j> A-d^^ ^ ^exp (^Ig J B dS^"^ , (107) 

where / is an integer and the integral is over the area E bounded by the loop 
C. We have 

Wi[C] = ^JJe''™"^ Z^^ J DAexp (^-AQ-^A + ilg J B ds"^ (108) 
The second factor is a Gaussian integral, which gives the factor 



Wa - exp ^ d^x d% v2G(x - y; 

In the limit of large E the leading behavior of the exponent is 



(109) 



;2 ;2 

g^j: lim k^G(k) = o^^E . (110) 

This gives the area law with the string tension^ 

a^^^g'm. (Ill) 

The first factor in Eq. (108) is different from unity only for odd I. It can be 
easily calculated to be 



+ 2z d^x -.COS (^xix) — a{x)j 



(112) 

'^As was shown in Ref. |13| the dependence on I in this formula is incorrect. The correct 
result is (T (X Z rather than a oc P. For multiply charged Wilson loops the nonlinearities of 
the compact theory are important and the Gaussian variational ansatz may be inadequate. 
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where a(x) vanishes for x outside the loop and is equal to tt for x inside the 
loop. At weak coupling we expand around a classical minimum of the exponent. 
Recall that the inverse propagator is nonzero at zero momentum Eq. (^0|). 
This dictates the leading order solution xi^) = and 



Ill 



, smce g m • 



(113) 



-7r/4/ 



This is a sub-leading correction to the string tension 
and z cx m? <C g^m. 

The behavior of the spatial Wilson loop suggests that the theory is con- 
fining with the string tension related in the expected way to the dynamically 
generated scale, a oc g^m. However, the spatial Wilson loop does not directly 
give the potential between external charges. Although in the Euclidean for- 
mulation there is no difference between spatial and time- like Wilson loops, 
in the Hamiltonian approach one should be more careful. In particular, the 
Hamiltonian formulation docs not preserve the Lorentz symmetry explicitly 
It is therefore important to calculate the potential between external charges 
directly. Our variational ansatz can be extended to calculate explicitly such a 
potentialE3 



5.6 Interaction potential between external charges 

In the charged sector the gauge invariance condition (|7^) is replaced by 



■(/"[A + VA] = exp {^g ^ p^X^ 



(114) 



where pn is a fixed, integer background charge distribution. To calculate the 
interaction potential between two charges we will later take it to be a well- 



separated dipole. A wave function that satisfies (114) is 



] = /M<^n]exp -i^4f™)G-i(r-s)4r) 



exp {^g ^ Pn(t>r}j 



The shifted field in ( |115| ) is defined along the lines of (pq), 
^(0,™) = ^ _ V(/, - • m . 



(115) 



(116) 
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but now it is necessary to define the vortex field in a singular gauge. Like 
A^, the field A^ satisfies 



fV X A'^) 



27r 

g 



-S hnJ 



bn',0 ■ 



(117) 



While A^ is divergenceless we take for A^ the solution where A is non-zero 
on the links dual to a string extending from (the plaquette of the vortex) to 
X — +00, 



a: 



ny 



0, 



^ for n.j: >0, ny = 0, 
otherwise. 



(118) 



Our reason for using A^' rather than A^ in (115) is one of locality. A dipole 
can be created in the vacuum ( |69| ) by the string operator exp [ig ^(z)] , 
which places sources at x and y. This operator creates a string of electric flux 
taking integer values along the string, the most local way of preserving Gauss' 
law when creating a dipole. A shift of A by A^ commutes with this string 
operator, but a shift by A^ does not. Shifting by A^ will create a non-local, 
transverse electric field with fractional flux in addition to the string. In this 
light it appears unavoidable that the introduction of dynamical charges will 
immediately lead to a nonlocal and non-Lorentz-invariant theory. For more on 
this point see Ref. |^ 

Another way of looking at it is to note that the difference between A^ and 
A^ can be absorbed by a shift in the integration variable 0r by -9{r — r')mr', 
where 9{r — r') is the angle that the vector r — r' makes with the x axis. Using 
A^ would therefore lead to an extra phase factor 



exp 



r,r' 



in the integral in (115). Under a shift of vortex density m this function is not 
invariant, but rather acquires a phase proportional to the charge density. The 
shift of vortex density can be viewed as a kind of large gauge transformation^ 
and therefore this wave function belongs to a sector of the Hilbert space with a 
position-dependent "0 angle." It is hard to imagine how this sector can define 
a local theory, especially in the presence of dynamical charges. 

Since A^ differs from A^ by a gauge transformation, it can be used in- 
terchangeably with it in the vacuum wave function (py). The distinction is 
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only meaningful in the charged sector, where it leads to a new vortex-charge 



interaction (see (123) b elow ) 
The wave function ( 



Hi 



satisfies the Gauss' law, but it turns out to be a 
very poor variational stateffij The reason is that it does not contain any extra 
variational parameter on top of the width of the Gaussian. Thus, the way 
the electric field, created by the external charges spreads in this state is not 
optimized. One can in fact calculate the distribution of the electric fieldtj 
This turns out to have a Coulomb-like profile. As a result, the energy of such 
a state is infrared divergent even for a finite size dipole. To do better we must 
introduce an extra variational parameter which could optimize the profile of 
the electric field. Thus, rather than using Eq. (115) we take as out state 



{'"„' } 



,{4>,m) 



exp ^ pn'/'n^ exp ^ ^ 



tin -f^n 



(119) 



We take the classical background field e to be transverse, V • e — 0, and we 
will treat it as an additional variational parameter in the charged sector. 
The normalization factor for the wave function (119) is 



ZAZ^[p]Zy[p] 



with 



Za = dot ttG 



z4p] = 



Zv[p] 



det47r^G 



1/2 



cxp -g'^p ■ -=7^0 ■ p 



(120) 

(121) 
(122) 



-i^E"V^(r'-s')r 



X exp ^ — i fer'^tr'^ ■ 

The new ingredient in ( |123| ) is a vortex-charge interaction, 
^Vmr' EE ^[pr6'(r-r') + '^^e^J^^{r~ T')ej{r) 



(123) 



(124) 
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The vortex interaction potential D is again given by Eq. (pO[). 

Let us now calculate the expectation value of the electric field in this state. 
Straightforward algebra gives 

(£;.>=g|^p-^G-i^(m)+e,. (125) 

The energy of this state is 

+^ j dPxSyd'^G-\x~y) {m{x)m{y)) 
-^Jd^x (Re(e-(^))-l) 

+ tCG-'{m) + lJd'xe\ (126) 
where the penultimate term contains the potential C{x), defined via 

e^ = ^^vdjC ■ (127) 

To calculate correlation functions of the vorticity m, we again introduce 
the dual field x as in Eq. (|8^) and obtain the relations 

(to) = -2ig^D-'{x), 

{m{x)m{y)) = 2g^ D'^ ~ ^g^ {D-\{x)D-\{y)) . (128) 

The duality transformation now results in the following Lagrangian for the 
dual field X: 

£./xi5-'x-2./.^x :»,(,(.)+<-..). ^ (129) 

The notation 9p represents the convolution of the source distribution p(v) with 

e{v-Y'). 

To first order in z we obtain 

{m{x)) = —2izsm{9p{x)~({x)), 
{■m{x)m{y)) — 2z cos{6p{x) — C{x))S'^ {x — y), 

{^e"'"^^")^! = exp[-4zcos(6i/9(a;)-C(a;))]. (130) 
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The e-dependent piece of the energy is 



AE 



2 2 

m g 



J (fx[l 



cos{9pix) - Cix))] 



+z / d^xd^yC{x)G~\x~y)smi9p{y)~ay)) 



(131) 



The quantity (131) is to be minimized with respect to ^. It is obvious imme- 
diately that at large distances from the sources C 9p, so that the energy is 
infrared finite. Noting that the second term is of order relative to the first, 
we drop it from now on. The minimization equation for C then becomes very 
simple, 

V^C - sin(C - 61/9) = . (132) 



To study Eq. (132) it is convenient to define C, — Q — Op which satisfies the 



sine-Gordon equation with a singular source term, 

V^C-m^ sinC = 5. 



(133) 



The source term S consists of a dipole layer along the line between the external 
point charges. For charges separated by a distance L, we have 



5(x) = 2^<5'(x2)77( ^ 



Xi 



Xi 



(134) 



wher e rj i s a step function. When L is much larger than 1 / gm the solution of 
Eq. (133) can be found in the region ~L/2 ^ a;i <C L/2. In this region C, is 
approximately independent of xi and satisfies the one-dimensional sine-Gordon 
equation in X2, 

^ m^sinC^S*. (135) 



dxl 



The singular source S creates a discontinuity of 27r at X2 = 0. The solution 
for a;2 > is hence half a sine-Gordon soliton, with — n dX X2 — 0+ and 
C, ^ Q as X2 ^ oo. For 0:2 < it is half the anti-soliton with Q = — tt at 
X2 = 0— and ^ — > as 2:2 — > ~oo. At distances greater than 1/m from the Xi 
axis, C, vanishes exponentially. Therefore, we have 9p exponentially. The 
vortex density m{x) according to Eq. (13C) vanishes exponentially. Referring 
to Eq. (12E) we see that the total electric field vanishes exponentially outside 
the region of width l/m, indeed forming a flux tube of thickness 1/m. This 
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is perfectly in accord with our expectation for the field profile in a confining 
theory. 



In Eq. (131) we see that the energy density of the flux tube is proportional 



to the energy of the sine-Gordon soliton solution with the proportionality co- 

2 

efficient The string tension is therefore 

2 

a ~ — g^m . (136) 

This is consistent with the result of the calculation of the spatial Wilson loop 
in the previous subsection. We have thus established that the electric field in 
the vacuum of compact QED is squeezed into flux tubes, and that the energy 
of a large dipole is proportional to its length. 

5. 7 Relation between the Lagrangian and the Hamiltonian pictures 

Some words should be said about the relation of the two-dimensional Hamilto- 
nian calculation to the three-dimensional Euclidean path integral. The vacuum 
wave functional of the theory can be represented in path integral formalism. 
To get the vacuum WF 'i'[A] one should calculate the path integral over the 
fields A{x,t), with t varying from — oo to 0, with the boundary condition 
A{x,t = 0) = A{x). To be more precise, in calculating the VEV of some 
operator 0{t = 0), one should split the time coordinate of the plane with the 
time coordinate of the operator, so that one considers 'i>[A{x,t ~ — e)] and 
'if*[A{x,t = e)] in the limit e — > 0. 

The basic objects that appear in the Euclidean path integrals are monopoles, 
which in 3D are not propagating particles, but rather instantons. When 
described in terms of the vector potential, or noncompact field strength, a 
monopole has a Dirac string attached to it. It is clear that the vortices (anti- 
vortices) in Eq. (|69| ) correspond precisely to the intersections of the Dirac 
strings of the 3D monopoles (anti-monopoles) with the equal time plane at 
t — 0. The positions of the Dirac strings are not physical in the compact the- 
ory, and only the position of the monopole itself is gauge invariant. In fact, for 
all monopoles that do not sit in the infinitesimally thin time slab between the 
planes t = — e and t = e, one can always choose the direction of the Dirac string 
such that it does not intersect the two planes. This precisely corresponds to 



the shift of the integration variable in going from Eq. (jTIj) to Eq. (72). The 
combination that enters this path integral nontrivially is the difference in vor- 
ticity between ^ and 5**. At points where both m (which corresponds to ^P) 
and to' (which corresponds to ^>*) are nonzero but equal, the difference van- 
ishes. This is the situation when a Dirac string intersects both planes t = ±e. 
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When a 3D monopole sits in the slab, only one of m or to' does not vanish 
and so their difference is nonzero. The summation over m in Eq. (^2|) can be 
interpreted therefore as the direct contribution to the partition function due 
to the monopoles at precisely the time t — 0. 

The fact that in this way one sees directly only the monopoles at t = 0, 
does not mean of course that other monopoles are not taken into account in 
this approximation. Indeed, the "bare" interaction potential between the t = 
monopoles is D{x). In the best variational state it is already short range, as 
follows from the solution for G{k) (see Eq. (|lOl| )). This is in accord with the 
3D picture, where the 3D monopole gas produces screening. Obviously, if one 
only looks at the thin slab, every monopole there will have an anti-monopole 
partner, which sits nearby (inside the screening length) in the third direction. 
The 2D monopole gas will therefore be screened by the 3D interaction, even 
before the interactions of the 2D monopoles between themselves is taken into 
account. This is perfectly consistent with our calculation. It is interesting to 
note, that even though this 2D interaction produces additional screening (the 
cosine term in the effective theory (pO|)), it is the 3D screening that is respon- 
sible for the area law of the Wilson loop, as is clear from the calculation of the 
string tension in Eq. (lOS). In fact, if one takes the noncompact expression 
G~^{k) = \k\, even the subleading term in the area law in the Wilson loop 
(113) disappearsEd 



6 The deconfining phase transition 

Now that we understand confinement in some detail, it is time to ask how it 
disappears. The deconfining phase transition in QCD at finite temperature 
has been the subject of numerous studies in recent years. Some aspects of the 
physics of the high temperature phase appear to be perturbative and can be 
studied in a controlled way at asymptotically high temperatures. However, in 
the phase transition region itself the QCD coupling is large, and the physics is 
dominated by the nonperturbative soft sector. Analytic understanding of this 
region is therefore extremely difficult. It is thus again useful to fall back on 
our toy model where a similar deconfining phase transition is within the weak 
coupling regime and can be studied analytically. 

In this section we discuss in detail the physics of the deconfining phase tran- 
sition. Analysis of the deconfining phase transition in the compact U(l) theory 
was performed in Ref. |2^. The SU(2) Georgi-Glashow model was studied in 
Ref. Interestingly enough, although at zero temperature the differences be- 
tween these two models are minute, the behavior close to the phase transition is 
very distinct. So much so that the dynamics that drives the phase transition as 
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well as the universal characteristics of the transition are completely different. 
In particular, whereas in the compact U(l) theory the transition is apparently 
driven by the binding of monopoles, in the Georgi-Glashow model it is due to 
the plasma effects of charged excitations {W bosons). In the lattice framework 
the monopole binding in the compact U(l) theory has been recently studied 
in Ref. |27 . One word of caution though, is that the compact U(l) theory does 
not have a proper continuum limit. Thus, it is not entirely clear whether the 
monopole binding mechanism is relevant to continuum field theory. In this 
review therefore we will concentrate on the Georgi-Glashow model. 



6.1 The monopole binding 

Since at zero temperature the monopole contributions are the only relevant 
ones for confinement, one may be tempted to assume that also at finite, low 
enough temperature all other effects are unimportant. The expected value of 
the transition temperature, as we shall see in a moment, is of order g^, which 
is much smaller than any mass scale in the theory except for the photon mass. 
One could therefore start with the working hypothesis that the single self- 
interacting photon field (or, equivalently, the monopole ensemble) should be a 
valid description of the phase transition region. This scenario was discussed in 
Ref. ||. 

The physics here is simple. The first thing to note is that at finite temper- 
ature the interaction between monopoles is logarithmic at large distances. The 
reason is that the finite temperature path integral is formulated with periodic 
boundary conditions in the Euclidean time direction. The field lines are there- 
fore prevented from crossing the boundary in this direction. The magnetic 
field lines emanating from a monopole have to bend close to the boundary and 
go parallel to it. So effectively the whole magnetic flux is squeezed into two 
dimensions. Qualitatively the situation is shown in Fig. 3. 



The length of the time direction is /? = l/T, and thus the field profile is 
clearly two-dimensional on distance scales larger than /?. Two monopoles sepa- 
rated by a distance larger than f3 therefore interact via a two-dimensional rather 
than a three-dimensional Coulomb potential, which is logarithmic at large dis- 
tances. Since the density of the monopoles is tiny cx ^, the monopole gas 
becomes two-dimensional already at extremely low temperatures T cx f^/^. 
The strength of the logarithmic interaction is easily calculated. The magnetic 
flux of the monopole should spread evenly in the compact direction once we are 
far enough from the monopole core. The field strength should have only com- 
ponents parallel to the spatial directions. Since the total flux of the monopole 
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l/T 



Figure 3: The field of a monopole-instanton at finite temperature. 

is ^ir/g, the field strength far from the core is 

F - 1^ 
g 

and thus the strengtiLof the infrared logarithmic interaction is j . 

It is well knowrO that the two-dimensional Coulomb gas undergoes the 
Berezinsky-Kostcrlitz-Thouless (BKT) phase transition. In the usual Coulomb 
gas, where the strength of the interaction A does not depend on temperature, 
the particles are bound in pairs at low temperature and unbind at high tem- 
perature T > Tbkt = SttA, where the entropy overcomes the energy. This is 
the standard BKT phase transitions which determines the universality class 
of U(l) symmetry restoration in 2 dimensions. In the present case the situ- 
ation is reversed, since the strength of the interaction itself depends on the 
temperature. At low temperature the interaction is weak, and therefore the 
particles (monopoles) are free. As the temperature grows, the interaction be- 
comes stronger until at Tbkt — 27rTgj^rp/g^ the energy overcomes the entropy. 
Above this temperature, the monopoles bind into neutral pairs. Thus, based 
on this simple picture one can expect the theory to undergo a BKT phase 
transition at 

Tbkt = |-. (137) 

ZTT 

Below this temperature the photon should be massive, while above this tem- 
perature it should be massless, since the cosine term in the Lagrangian (^) is 
irrelevant. 

This is precisely the logic followed in Ref. One can go further and 
perform more quantitative calculations using the dimensionally reduced ver- 
sion of Polyakov's effective Lagrangian (|l4|) . Dimensional reduction should be 
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perfectly valid in this theory since the critical temperature is much larger than 
the photon mass. The effect of the monopoles at finite temperature is thus 
contained in the two-dimensional model 

This theory, in full agreement with the previous discussion has a BKT phase 



transition at the temperature given by Eq. (137), above which the "dual pho- 
ton" field (j) becomes massless. 

Moreover, since the two-dimensional sine-Gordon theory is exactly solv- 
able, one can calculate various correlation functions. In particular, the string 
tension was calculated in Ref. and was found to exhibits quite strange 
behavior. At low temperatures it follows an expected pattern — that is, it 
decreases with temperature. It becomes extremely small at temperature much 
smaller than Tbkt; however, just before Tbkt it starts rising sharply and at 
Tbkt actually diverges. 

Although the logic leading to the previous discussion seems quite sound, 
there are several puzzles that arise regarding the results. 

• The photon becomes massless in the high temperature phase. In other 
words, the correlation length in some physical (gauge invariant) channel be- 
comes infinite at high temperature. This contradicts our physical intuition, 
since one expects that at high temperature any physical system described by 
a local field theory should become maximally disordered with zero correlation 
length. 

• The phase transition is of BKT type, and therefore belongs to the U(l) 
universality class. On the other hand, we knowO that the global symmetry 
that is restored at the phase transition is the magnetic Z2, and we expect 
the universality class to^be 2D Ising. This is also expected from the classic 
universality argumentstJ 

• The divergence of the string tension at the critical point looks unphysical. 

• Finally, as we explained in the previous sections, the effective theory (14) 



does not allow charged states. It is then difficult to understand in what sense 
the high temperature phase can be viewed as deconfined. 



6.2 Magnetic symmetry restoration and the charged plasma 

Another line of argument starts with the realization that confinement is tan- 
tamount to the spontaneous breaking of the Z2 magnetic symmetry. Indeed, 
let us now take Eq. (^ ) as our starting point. What does one expect from 
the phase transition in a simple scalar theory of this type? Firstly, clearly one 
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expects the eritieal temperature to be of order of the expectation value of the 
scalar field, and therefore parametrically of order g^. At finite temperature 
one expects the generation of a positive thermal mass proportional to the tem- 
perature and to the coupling constant. Thus, the thermal contribution to the 
effective potential should be of the form 

5tC = x\TV*V . (139) 

Ax T — /Att^x the total mass term becomes positive, the VEV of V van- 
ishes and the phase transition occurs. One can estimate this temperature 
more precisely using the following simple argument. Let us neglect for now 
the monopole induced term. Then we are dealing with the XY model at finite 
temperature. The dimensional reduction should again be a valid approxima- 
tion, and thus essentially we have to analyze a 2D model. Now a 2D XY model 
can be mapped into a sine-Gordon theory of a dual field Performing this 
dual transformation^, we find the Lagrangian 

c^^{^^xf + ^l cos x- (140) 

where /Lt is the fugacity of the vortices in the XY model. This sine-Gordon 
theory has a phase transition at 

T.y ^ ^. (141) 

Therefore, we may expect the magnetic symmetry restoring phase transition 
at TxY- Although parametrically this temperature is of the same order as the 
BKT transition temperature discussed in the previous subsection, it is four 
times lower. Another important difference is that the nature of the phase 
transition is in fact completely different. The phase transition in the model 



(140) is due to the unbinding of vortices of the field V. Above Txy the vortices 
are in the plasma phase. These vortices are, as discussed earlier, precisely the 
charged bosons of the original Georgi-Glashow model. Thus, this phase 
transition is just what one would naturally call the deconfining phase transition 
and the vacuum above Txy is the charged plasma. 

Note that in this discussion we have neglected completely the effect of 
monopoles, that is the last term in the Lagrangian (|2^). Interestingly enough 
we see that the plasma phase is reached at a temperature which is much lower 
than Tbkt discussed earlier, and thus the monopole binding is irrelevant for 

"^To fix the normalization of the kinetic term we should bear in mind that a vortex corre- 
sponds to a 2iT jump in the field x = 2(/> rather than the field </> in Eq. (fl 
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the dynamics of the deconfinement. This is not to say that the presence of the 
monopoles is irrelevant altogether. Clearly, omitting the monopole induced 
term, we enlarged the symmetry of the system from Zi to U(l). Hence the 
effective XY model description and the U(l) universality class predicted by 



Eq. ( 140 ). The analysis as it stands is correct for noncompact electrodynamics 
with charged matter, but not for the Georgi-Glashow model. Next we remedy 
this problem. 

The effect of the monopoles can be qualitatively understood in a simple 
way. Their presence leads to a confining potential between the charges, that 
is a linear interaction between the XY model vortices. Thus, the monopoles 
suppress the variation of the phase of the field V except inside the confining 
string, the width of which is given by the inverse photon mass d — 1/M cx 
Whether the presence of the monopoles is important at the would be phase 
transition point Txy depends crucially on the density of charged particles. If 
the density of the charged particles is very high, so that the average distance 
between them is smaller than rf, the presence of the monopoles is immaterial 
since the phase of V is disordered already on the short distance scale. However, 
if the density of charges is low, so that the distance between them is larger than 
d, the presence of the monopoles will suppress the phase transition. In this 
case at Txy the vacuum will not be disordered, but will look like a dilute gas 
of charged particles with strings between them. The actual phase transition 
will then occur at a higher temperature, where the average distance between 
the charges equals the inverse photon mass. 

Interestingly enough in the present case this temperature turns out to be 
twice the value of Txy and therefore still much lower than Tbkt- We will 
show this rigorously in the next subsection. Qualitatively though this is easy 
to understand. The density of charges is proportional to their fugacity 

fi(xe^i^. (142) 

This should be compared to the square of the photon mass, which is given by 
the zero temperature monopole fugacity. In a theory with very heavy Higgs 



C oc e . (143) 

The two are equal (up to subleading corrections) at 



^2 



Tgg = f^. (144) 
More generally, when the Higgs is not infinitely heavy, the monopole fugacity 
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is smaller and thus the transition temperature will be lower, 




(145) 



For a lighter Higgs, Tc gets closer to Txy, but is always greater, so that 
TxY <Tc < Tbkt- 

In the next subsection we will present a more complete analysis based on 
the renormalization group and exact bosonization. This analysis confirms the 
simple picture presented here. Thus, all the puzzles raised in the previous 
subsection simply disappear: 

• The photon never becomes massless. Even without the monopoles in 
the plasma phase it acquires the Debye mass given by the cosine term in Eq. 
( |140| ) . This mass rises with temperature, and thus the physical correlation 
length decreases. 

• Since the monopole term is still relevant at the phase transition, the 
universality class must be Z2- We will see this explicitly in the next section. 

• The analysis of the previous subsection is only valid below Txy, and 
thus the divergence of the string tension at Tbkt > Txv has nothing to do 
with the physics of the Georgi-Glashow model. 

• Finally, the phase transition is driven primarily by the unbinding of 
charged particles and thus indeed has all the flavor of a deconfinement transi- 
tion. 

Let us now turn to a more quantitative analysis of what is stated above. 
6.3 The renormalization group analysis 

For a more formal analysis of the phase transition we find it convenient to use 
the sine-Gordon formulation in terms of the phase field x modified to take into 
account explicitly the finite probability of the appearance of vortices. This 
Lagrangian has the form 



where ( is related to the monopole fugacity C, — ^/T and x is the field dual to 



C = 



87r2 




{dt,xf + C,cos2x + ^i cos X , 



(146) 



X, 



■id^X 




(147) 
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The way to derive this Lagrangian is as follows. The partition function of the 
sine-Gordon model in the presence of one vortex is 



Z{x) = j Z?[x]exp ~ j d 



'^y^2f (^mX - Jm (y. ^)f + C cos 2x 



The "external current" is 



(148) 



(149) 



with C a curve that starts at the location of the vortex (the point x), and 
goes to infinity, and n^^ is the unit normal to this curve. The insertion of this 
current forces the derivative of x to have a discontinuity across the curve C, so 
that X jumps by 27r. This forces x to have one unit of vorticity concentrated 
at the point x. Note that even though explicitly depends on the curve C, 
the partition function itself does not, since changing the integration variable 



X{x) x{x) + 2tt, 



X E S , 



(150) 



where the boundary of 5 is C — C", is equivalent to changing C into C" i n the 
definition of the current. The extra linear term in the exponential in Eq. (148) 
is 



IX 



9 

2ttT 



/ dx^e^^d^x, 
Jc 



(151) 



which is equivalent to Eq. ( 147 ). An anti- vortex at y is obviously created by 
—jfj,. To create several vortices one just inserts an external current which is 
the sum of the currents which create individual vortices. 

A dilute ensemble of vortices and anti-vortices with (small) fugacity /i is 
then given by 



^ = / n,dx,Il,dy,Z{x,,y,). 



(152) 



The summation over the number of vortices and anti-vortices can be easily 
performed leading to the partition function with the Lagrangian (146). The 
constant fi is the vortex fugacity scaled by the effective ultraviolet cutoff im- 
posed on the integration over the coordinates. The vortex fugacity of course 
is none other but the fugacity of the charged W, 



(153) 



The cutoff a is related to the Compton wavelength of the W boson, but a 
more careful determination of it should take into account the fact that in the 
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process of dimensional reduction, all modes with frequencies above T have 
been integrated out. We will not attempt the determination of a, but only 
note that it is some combination of the scales Mw and T oc g^, and as such 
its value always plays a role secondary to the exponential factor of fugacity. 

An alternative way to derive the Lagrangian ( |l46|) is to start directly from 
the effective theory ( p8|) with the extra Skyrme term (^5|). In the nonlinear 
cr-model limit one can cleanly separate the phase of the field V into a smooth 
part X Etnd the vortex contribution. The Skyrme term then is proportional to 
the energy of one vortex and, with the Skyrme coupling chosen the way we 
discussed in Sec. 2, is just equal to Mw for a one vortex configu ratio n. The 
dilute vortex gas approximation then reproduces the Lagrangian ( |146[ ) . 

Since both ^ and fj, are small, the importance of different terms in the 



Lagrangian (146) is determined by their respective conformal dimensions cal- 
culated in the free theory. The total conformal dimensions of the operators 
: cos2x : and : cosx respectively, are 

A, = — and A, = ^. (154) 

Thus, at low temperature the charges are irrelevant (A^ > 2), the monopoles 
are relevant (A^ < 2) and the theory reduces to the sine-Gordon model. At 
the temperature T = (?^/47r the two operators become equally relevant, since 
their conformal dimensions are equal, 

A^ = A^ = 1. (155) 

The phase transition point can be determined from the structure of the 
fixed points of the renormalization group equations. The critical point is the 
infrared unstable fixed point of the renormalization group flow. The renor- 
malization group equations of the sine-Gordon theory have been studied per- 
turbatively in the literature both in the absencdj and in the presenceEJ of the 
vortices. In terms of the dimensionless parameters 

t^^T, jl:^fia^, C = Ca^ (156) 
9 

the lowest order renormalization group equations read 

^ = n'ip'-t^C'l (157) 

fx - (2-0 C. (159) 
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The fixed point structure of these equations is simple. 

1. The point Tq 

i = 0, /i = 0, C = oo (160) 

is clearly the zero temperature fixed point. Here the long distance physics is 
dominated by the monopole induced mass term. 

2. The point T<x, 

t = oo, jl = oo, C = (161) 

is the high temperature fixed point. Here the infrared properties are deter- 
mined by the charged plasma eflfects. 

3. The point Tqc 

t = l, il = C, /i = oo (162) 

is the infrared unstable fixed point. This is precisely the critical point that cor- 
responds to the deconfining phase transition. At this point both the monopole 
and the charge plasma induced terms are equally important. 

Note that the fixed line t > 2, C, = oi massless theories which corre- 
sponds to the "confined monopole plasma" is not present, since for i > 2 the 
charged plasma induced mass term is strongly relevant. The same is true for 
the would be fixed line t < 1/2, p. = 0, which is present in the absence of the 
monopoles and describes the low temperature massless phase of noncompact 
electrodynamics . 

It is instructive to see how the renormalization group equations formalize 
our qualitative arguments of the previous section. In particular, they make 
clear the role of the point t = 1/2 which, in the absence of monopoles, would 
be the point where the charged induced term becomes relevant. If one starts the 
evolution from the initial condition t = 1/2, /i ^ C, the running temperature 
will increase, and jl will grow into the infrared, while ( will grow for a while until 
the running temperature reaches two. From this point on, ( will decrease and 
the system will approach the high temperature fixed point. This corresponds 
to the situation where at the would be critical point t = 1/2, the density 
of charged plasma is so high that the mean distance between the particles is 
smaller than the width of the confining string. 

On the other hand, if the initial condition is /i < i^. the temperature 
starts decreasing, making the coupling jl immediately irrelevant. Thus, jl 
monotonically decreases to zero and does not significantly affect the flow of 
the other couplings, which steadily flow to the zero temperature fixed point. 
This holds in the situation where at t = 1/2 the density of charged plasma is 
low, which is indeed true in our model. 
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Interestingly enough, at < = 1 the initial conditions in our model in the 
case of a very heavy Higgs are such that the renormalization group flow starts 
almost exactly in the region of attraction of the fixed point Tqq. At this 
value of temperature the monopole and charge fugacities are equal. The only 
difference between jl and C is then in the prefactors, which are possibly different 
combinations of Mw and g^. If /i and ^ were exactly equal, the only thing that 
happens along the flow is that their values grow, but t does not change. Due 
to the small difference their values are equal at a temperature which slightly 
differs from 1, 



T 



9^ 
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l + O 



\Mw) 



(163) 



Thus, the initial temperature for which the system is in the region of attraction 
of Tgg is slightly different. We conclude that the critical temperature of the 

2 

Georgi-Glashow model is indeed given by Tgg = f- , up to corrections of order 



w 



6.4 The universality class of the phase transition 

Since the deconfining phase transition is due to the restoration of Z2 symmetry, 
we expect that it is in the two-dimensional Ising universality class. This can 
be shown exactly by studying the theory ( 146| ) at the fixed point t = 1. The 
following discussion follows closely Rcf. [See also Chapter 21 in the book 

The theory described by Eq. ( |l46[ ) can be fermionized by using the stan- 
dard bosonization/fermionization techniques. Since at i = 1 both cosine terms 
have dimension 1, the resulting fermionic theory is a theory of free massive 
fermions. In our notations, the Dirac fermionic field is defined as 



V'fl — a ^^^i exp 



^+2 



tpL — o, exp 



(164) 



The kinetic term in Eq. ( |146| ) then becomes the kinetic term of the field tp, 
while the cosine terms become 



H.c. 



a ^ cos 2x^ i[>pli'>pL 

a^"^ cosx = ~ H.c.]. (165) 

Thus, the mass term in the fermionized Lagrangian is diagonalized by intro- 
ducing the real Majorana fermions 

ip + tp^ -0 — -0'!' 



V2 



iV2 



(166) 
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The mass of the fermion p is /ia + Ca, while the mass of a is fia — (a. For /i = C 
the model contains one massive and one massless fermion. At the fixed point 
// — > oo the massive fermion decouples. Thus, at the point Tqg the theory 
is that of one massless Majorana fermion. It is well know that this theory 
precisely describes the critical point of the 2D Ising model. 

Thus, we indeed see that the phase transition is in the universality class 
of the 2D Ising model. 



7 Vortices and monopoles at high temperature 

It is interesting to see how the behavior of magnetic vortices and the monopoles 
is affected by the deconfining phase transition. We will do this with the help 
of the effective action written in terms of the zeroth component of the vector 
potential Aq. The relation between our analysis, which was performed for the 
phase of the vortex operator x and this standard procedure is not difficult 



to see. The free energy of the charged bosons in the Lagrangian (146) 
is given by the term cosx. On the other hand, this same free energy in the 
standard calculation is represented by the insertion of the Polyakov line with 
charge two. Thus, we identify the Polyakov line with the exponential of the 
dual field 

P = exp (^^x) , (167) 
and the dual field x with the Abelian vector potential 



(168) 



With the identification in 



Up to the monopole induced term, the Lagrangian (146) is equivalent to the 
sine-Gordon theory of the field x in Eq. 
Eq. ( |168| ), this is 

^-Xd^Aof 



H cos 



25 

T' 



Ac 



(169) 



The monopole induced term also can be written in terms of the vector potential. 
Its form in terms of the Polyakov line P is precisely the same as the form of 
the charge plasma induced term in terms of the vortex field V. Recall that 
the origin of the plasma induced term cosx in Eq. (146) is the dimensionally 
reduced "Skyrme term" of Eq. (^5|). Thus, just like we derived the Sine-Gordon 
Lagrangian (146) starting from the effective Lagrangian for the vortex field, 
we can follow the same steps backwards, but this time expressing everything 
in terms of P. Due to this duality between V and P we conclude that the 



48 



monopole induced term is the dimensionally reduced Skyrme term for the 
Polyakov line. The fuU effective Lagrangian for the vector potential is then 



£ = j;{d,Aor + Mcos ( ^Ao ) + £^ In ( C) {e,,d,Pd,P 



T 



47r2 



(170) 



The monopole in this Lagrangian is represented by a vortex of the field 
with unit vorticity. The coefficient of the Skyrme term in Eq. (170) is such 
that the action of such a unit vortex is equal to the action of the "core" of the 
monopole. 

The cosine term in this expression is the potential for induced by the 
nonvanishing density of the charged particles in the thermal ensemble. Natu- 
rally, it contains the Debye "electric" mass term for Ag and also higher inter- 
actions. Note that, as opposed to strongly interacting theories, where similar 
effective Lagrangians have been derived only in the derivative expansion, the 
Lagrangian (170) is valid on all distance scales longer than 1/T. Thus, in prin- 
ciple, it can be used to calculate correlation functions in a large momentum 
range. However, if one is interested in the long distance behavior of the cor- 
relators, at temperatures above TqGj the higher derivative Skyrme term can 
be neglected. This is in accordance with our analysis of the previous section, 
which showed that the monopole induced term is irrelevant above Tqg- 

Consider now the correlation function of two vortex operators. Recall 
that in the confining phase the vortex operator has a nonvanishing expecta- 
tion value since the magnetic Z2 symmetry is spontaneously broken. In the 
deconfined phase we expect the correlation function of two vortex fields to 
decay exponentially at large distance. 

Consistent with the magnetic symmetry breaking, at zero temperature we 
have 



{V{x)V*{y)) = ^e^V 



\{x{x)x{y)) 



Stt- 



■ exp 



167r^ 
7^ 



Y:i{x - v) 



(171) 

where ¥3(2; — y) is the 3d Yukawa potential with the mass M . At temperatures 
below the phase transition (T < Tqq) but high enough so that the dimensional 
reduction is valid the infrared asymptotics (at distances |a; — j/| > 1/T) is 
instead 



(^(^)^*(y))-|^exp 



167r^T 



Y2{x~y) 



(172) 



with Y2{x — y) the two-dimensional Yukawa potential with the mass which 
includes the effects of integrating out the nonzero Matsubara modes, as calcu- 
lated in Ref. |5|. This expression follows from Eq. (146) neglecting the cosx 
term, which is indeed negligible in the infrared. 
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At high temperatures T > Tqq it is convenient to use Eq. (170) with the 
omission of the third term. The calculation of the correlator of the vortex 
operators proceeds along the lines of Refs. 0. The insertion of V{x) 
and V*{y) creates the Z2 domain wall stretching between the points x and 
y. In terms of the sine-Gordon theory, Eq. ( |17G| ) this domain wall is just the 
kink, and thus the domain wall tension is equal to the soliton mass Mg. The 
correlator is then 

{V{x)V*{y)) = exp{-M,|a; - y\} (173) 

with the soliton mass 



M, 



-1 2r(p/2) 



2/i 



(174) 



where 



P 



9 



SttT - 52 ■ 



(175) 



Thus, we find that the correlation function of the vortex operators de- 
creases exponentially in the high temperature phase as it should in the phase 
with restored symmetry. Note that, were we to neglect the Debye mass term, 
we would have found that below Tbkt the correlator tends to a constant at 
infinity, while at T > Tbkt it decays at large distances, but only as a power, 
since the mass of the soliton vanishes. 

What about the monopoles? The discussion in the beginning of this section 
would suggest that they are "confi ned" logarithmically. However, it is easy to 
see from the effective action Eq. (17C) that the monopoles are bound by a 
(screened) linear potential. Consider a configuration with unit winding of the 
field P. Due to the potential term P^, the minimal action configuration cannot 
be a rotationally symmetric hedgehog. Such a configuration would "cost" 
action proportional to the volume, since the field P would be away from its 
vacuum value everywhere in space. The best one can do is to have a quasi one 
dimensional strip in which the winding is concentrated, while everywhere else 
in space P would be equal to 1 (or —1). This configuration is schematically 
depicted in Fig. 4. 

The width of the "dual confining string" must clearly be of order of Mjj^, 
while the action per unit length, s cx ^Md- Thus, the action of a sin- 
gle monopole diverges linearly with the size of the system. Obviously the 
monopole-anti-monopole pair separated by a distance L have action Ls. When 
the distance is large enough, another pair can be produced from the thermal 
ensemble, screening the linear potential. The critical distance at which this 
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Figure 4: The string-like configuration of the Polyakov line P that accompanies the 
monopole-instanton in the high temperature phase. 



happens is determined by 



ic. = 21ogC = ^, (176) 

or 

L. oc ^M-. (177) 

Thus, as long as the temperature is much lower than Mw, the length of the 
"dual string" is much greater than its thickness. The "potential" between the 
monopoles is therefore linear, but screened at large distances, much the same 
as the confining potential in gauge theories with heavy fundamental charges. 
Thus, the behavior of the monopoles is in many senses "dual" to that of electric 
charges]^. 

Even though the monopoles are "confined" , their effects do not disappear 
in the hot phase. They, for example, are the main contributors to the correla- 
tion functions of the vortex operators, of the type V{x)V{y). This correlator 
would be vanishing if the magnetic symmetry were U(l) rather than Z2. It 
was shown, however, in Ref. p4, that this correlation function does not vanish 



■^Of course one always has to keep in mind that physically monopoles and charges are very 
difi^erent objects in this model; the charges are particles, while the monopoles are instantons. 
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due to the monopole contributions. In fact, the "diagonal" correlation func- 
tions are those of the parity even and parity odd operators = V + V* and 
V- =V-V*: 

{V+{x)V+{y)) = {aMo)^ eM-\x - y\^-} 
{V-{x)V-iy)) = (aMo)^ e^p{~\x ^ y\a+} (178) 

with 

a± = Ms ± CiaMD)"^'^ — - (179) 

na 

Thus, the correlation lengths in the scalar and pseudoscalar channels are not 
equal due to the monopole effects. At high temperatures this difference be- 
comes small, since C. <^ aMD <C 1, but it nevertheless remains strictly finite at 
any finite temperature. 

8 Hot confining strings: Deconfinement versus Hagedorn transi- 
tions 

In this section we discuss the deconfining phase transition from the perspective 
of the confining string. 

What happens to strings at very high temperatures is one of the big ques- 
tions of string theory. In the early days of dual resonance models and hadronic 
bootstrap, the exponential growth in the one-particle density of states 

p{m) - e*"" (180) 

led to the famous Hagedorn transition^ Tliia type of spectrum first arose in 
the context of statistical bootstrap modelsEjOEll and, for hadrons, such he^ 
havior indicates that they are composed of more fundamental constituents^ 
In fundamental string theories one finds the same kind of spectrum (see, for 
example, books |3^ ), and a search for hints to the existence of "string con- 
stituents is of great interest. What lies beyond the Hagedorn temperature? 
Is this temperature limiting or is there a high-temperature phase which re- 
veals the fundamental degrees of freedom? And is it true that for all types of 
string theories there is the same universal physics or is it that different classes 
may have totally different high-temperature behavior? A lot of effort has 
been invested into study of the Hagedorn transition in critical (super)strings. 
There is an enormous volume of literature on this subject, some references 
(but by no means all) can be found, for example, in Ref. For weakly 

coupled critical (super)strings the Hagedorn transition can be described as a 
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Berezinsky-Kostcrlitz-ThoulessH transition on a world shcctC^I'oHcJ It is due 
to the "condensation" of the world sheet vortices. It has been also suggested 
that the transition in some cases may actually be first order. Above the Hage- 
dorn temperature, the vortices populate the world sheet and we have a new 
phase. From the target space point of view we are talking about tachyonic in- 
stabilities for non zero winding modes in the imaginary time direction (thermal 
winding modes). 

In the more general case of interacting strings, one does not know with 
certainty what the fundamental degrees of freedom are and thus what their 
role is at high temperature. Recently it was suggested, in the framework of 
the Matrix Model description o£±he Hagedorn transition, that the fundamental 
string decays into DO Branesc3 Thus, it could be that DO branes are the 
fundamental degrees of freedom in the hot phase. 

In short, there are many interesting open string theory questions pertaining 
to the hot string. Since our simple 2+1 dimensional model is well understood, 
it is worth seeing what it has to say on these issues. 



8.1 The monopole binding as the Hagedorn transition 

As we saw in the previous sections, the dynamics of the deconfining phase 
transition is quite nontrivial, involving interplay between the charged degrees 
of freedom and the monopoles. To disentangle these effects in string language, 
we start our discussion by completely disregarding the charged particles. From 
the string point of view, this means that we neglect possible contributions of 
the heavy DO branes, and are thus entirely within the theory of closed strings. 
Naively one expects in such a theory the existence of a Hagedorn temperature, 
beyond which the string cannot exist. In an almost free string this temperature 
is of order of the string scale. One can visualize this phenomenon in simple 
terms. Consider a closed string of a given fixed length L. Let us calculate the 
free energy of such a string. The energy of the string is 

E = aL (181) 

The number of states for a closed string of length L scales exponentially with 
L 

N{L) ^ exp{aL}. (182) 

The dimensional constant a is determined by the physical thickness of the 
string. Imagine that the string can take only positions allowed on a lattice 
with the lattice spacing a. Then clearly the number of possible states is z"L, 
where z is a number of order unity, equal to the number of nearest neighbors 
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on the lattice.^ In this simple situation a = alnz. The only natural lattice 
spacing for such a discretization is the thickness of the string. For an almost 
free string the thickness is naturally the same as the scale associated with the 
string tension. Thus, the entropy is 



S = xy/^L (183) 
with X a number of order unity. The free energy is then 



F[L] = aL- xTy/^L (184) 

At the temperature 

Th = -V^ (185) 

X 

the free energy becomes negative, which means that strings of arbitrary length 
appear in the thermal ensemble in a completely unsuppressed way. The ther- 
mal vacuum becomes a "soup" of arbitrarily long strings. Thus, effectively, the 
"temperature dependent" string tension vanishes and it is not possible to talk 
about strings any more in the hot phase. For more details about this "random 
walk" description of hot strings, see Ref. ^ and references therein. 

The situation is very similar to the BKT phase transition, where the free 
energy of a vortex becomes negative at the critical temperature, and the vor- 
tices populate the vacuum in the hot phase. 

In the string partition function language this is just a restatement of the 
well known fact that the partition function diverges in a sector with the topol- 
ogy of a torus which winds around the compact Euclidean time direction. 
Fixing the unit winding in the Euclidean time physically corresponds to the 
calculation of the free energy in the sector with one closed string. The inte- 
gration over all possible lengths of the string is the cause of the divergence of 
the partition function at high temperature. 

The same physical effect exists for the confining Georgi- Glashow string. 
There is one difference, which however turns out to be crucial for the nature of 
the phase transition. This is that the thickness of the Georgi- Glashow string is 
not given by the string tension. Rather it is equal to the inverse photon mass, 
and thus 

a cx cx . (186) 

We disregard in this argument the fact that the string has to close in on itself. This extra 
condition would lead to a pre-factor with power dependence on L. Such a pre-factor is not 
essential for our argument, and we therefore do not worry about it. 
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We thus have for the confining Georgi-Glashow string 

Tff = - cx g". (187) 

a 

This is indeed the correct magnitude of the critical temperature as discussed 
in previous sections. 

The noteworthy feature of this formula is that the Hagcdorn temperature 
of the confining string is much higher than the string scale. This is easy to 
understand because the entropy of the thick string is much smaller than that 
of the free string due to the fact that high momentum modes of the confining 
string do not contribute to the entropy at all, as discussed earlier. Thus, one 
needs to heat the string to a much higher temperature for entropy effects to 
become important. 

Since we have completely neglected the possibility of the appearance of 
DO branes (charged particles), the transition we have been discussing is the 
string representation of the monopole binding transition in the Georgi-Glashow 
model. At the point where the monopoles bind, the photon of the Georgi- 
Glashow model becomes massless and thus the string tension disappears and 
the thickness of the string diverges. In the Hagedorn picture this is just the 
dual statement that the ensemble is dominated by infinitely long strings. The 
BKT nature of both transitions underlines this point. 

8.2 Vortices on the world sheet — open strings and charged particles 

Sometimes the Hagcdorn transition is discussed in terms of the vortices on the 
world sheet. What is the physical nature of these objects? 

Consider a string world sheet with the topology of a sphere with a vortex- 
anti-vortex pair, Fig. 5. When going around the vortex location on the world 
sheet, the compact coordinate xq varies from to /3. Since this is true for 
any contour of arbitrarily small radius, which encircles the vortex, this means 
that physically the location of the vortex in fact corresponds in the target 
space not to a point but rather to a line which winds around the compact 
direction. The vortex-anti-vortcx pair on the world sheet thus rc;prescnts an 
open string which winds around the compact direction. Figure 5 illustrates 
how the open string world sheet, equivalent to a cylinder is transformed by a 
conformal transformation into a sphere with two singular points — the vortex- 
anti- vortex pair. 

If the string theory in question does have an open string sector, then the 
configurations with arbitrary number of vortcx-anti-vortex pairs contribute to 
the finite temperature partition function. The Hagcdorn phase transition can 
then be discussed in this sector rather than in the closed string sector. Not 
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Figure 5: Representation of a vortex and an anti-vortex pair with a string stretched between 
them and the corresponding string world sheets. 



56 



surprisingly, the discussion is exactly the same as in the previous subsection. 

The fact that wc now have an open rather than closed string docs not change 
the entropy versus energy argument. One finds that at T < Th the vortices 
on the world sheet are bound in pairs. The corresponding target space picture 
is that the "ends" of the open strings are bound by a linear potential and 
therefore long open strings do not contribute to the thermal ensemble. In 
exactly the same way, long closed strings are also absent from the ensemble. 
At T > Th the vortices unbind and appear in the ensemble as a Coulomb gas. 
Thus, a typical configuration contains lots of open strings (as well as lots of 
arbitrarily long closed strings) since the energy of such strings is overwhelmed 
by the entropy. 

Note, however, that the existence of the transition in this context is entirely 
independent of the presence or absence of vortices. As we saw above, the 
transition can be understood purely on the level of the closed string. It is 
driven by the string fluctuations. Thus, even if the theory does not have an 
open string sector, the transition is still there. For example there are no open 
strings in the Georgi-Glashow model. Nevertheless, if we neglect the effects of 

the Hagedorn transition is still there and it coincides with the "monopole 
binding" transition. 

The actual phase transition in the Georgi-Glashow model is however not 
driven by monopole binding. Nevertheless, in the string language it is also due 
to proliferation of vortices on the world sheet. These are, however, vortices of 
a somewhat different type. The same conformal transformation that turned 
the open string boundary into a point can be used to turn the world line 
of a heavy particle into a point — let's call it brane. As discussed above, 
this is a (fundamentally) charged particle which couples to the Georgi-Glashow 
confining string. Thus, a string world sheet representation of a pair of particles 
with opposite charge is also a vortex-anti-vortex pair on a sphere. 

There are some important differences between these vortices and the ones 
that represent the ends of the open string. First, for an open string the non- 
compact "spatial" coordinates satisfy Neumann boundary conditions. Thus, 
even as winds, the other coordinates can take arbitrary values close to 
the vortex. One can see this from the action 

S = a j d'^^dc,x°dax° +a j d^^dax'dax' (188) 

where the dynamics of is absolutely unrelated to the dynamics of x^ . Thus, 
even if a vortex in sector can be considered as a boundary in a target space, 
there are free boundary conditions for the components x^. In other words, wc 
do not have any boundary action induced by vortices in a theory of closed 
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strings. 

For a very heavy brane the boundary conditions are of the Dirichlet type. 
Thus, Xi are constant close enough to the vortex location. Another difference is 
that, since the brane is an independent degree of freedom, in principle its mass 
is a free parameter. Thus, the fugacity of the brane vortex is an independent 
parameter, and physics may in principle depend on it. Of course, in the case 
when we have a finite mass for branes there is a nontrivial boundary action 
describing the massive particle. In that case one has pan-conformal boundary 
conditions compatible with a finite mass of a braneeZl 

The Georgi-Glashow confining string has neither open string sectors nor 
branes which are the sources of a single string. The dynamical objects which 
couple to the string are the charged particles, which have an adjoint charge 
and therefore are sources of a pair of strings. 

Each dynamical brane in the Georgi Glashow model has two strings 
emanating from it. Thus, a pair of branes propagating in compact imaginary 
time is conformally equivalent to a pair of "double vortices" as in Fig. 6. 

The singular points are still vortices as before, since going around such 
a point one travels once around the compact time direction. But now two 
string world sheets are permanently glued together at the location of a vortex. 
The fact that two world sheets are glued at the location of the vortex is the 
manifestation of the Z2 magnetic symmetry of the Georgi- Glashow model. 
The region of space between the two world sheets is separated from the region 
of space on the outside, reflecting the fact that it constitutes a domain of a 
different vacuum. 

At low temperatures such configurations in the thermal ensemble are rare. 
But at the critical temperature their density becomes so large that the ensemble 
is dominated by multi-vortex configurations. In the string picture these are 
configurations with multiple points of "gluing" . Importantly, the branes are 
not just pairwise connected by two strings, but, rather, form a network where 
all of them are connected to each other. Clearly in terms of entropy such 
configurations are much more favorable, and there is also no loss of energy 
when the distance between the branes is of the order of the string thickness. 

This string-based description of the transition has to be taken with a grain 
of salt. The relevant physics takes place on short distance scales — of the order 
of the thickness of the string. On these scales, as we explained above, the string 
modes are practically absent. Thus, even though we have showed the string 
segments on Fig. 6, these segments are so short that there is no string tension 
associated with them. Thus, the mechanism of the transition is essentially 
field-theoretical rather than string-theoretical. 

To understand this point better consider in a little more detail the thermal 
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Figure 6: See the text 
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ensemble of W^. The crucial point is that at distance scales d ^ A/^^ the in- 
teraction between them is Coulomb rather than linear. The gas of charges with 
Coulomb interaction has itself a transition into a plasma phase. This transition 
has nothing to do with the long range linear interaction and occurs at theJ^em- 
perature T^c which is four times smaller than the Hagedorn temperatureEj At 
this temperature the become "free" in the sense that they cease to care 
about the Coulomb part of the potential. The crucial question however is: 
how large is the density of at this point? If the density at T^c were large 
and the average distance between W's were smaller than M~^, the transition 
would actually occur at Tnc, since the long range linear part of the potential 
would be entirely irrelevant. As it happens, in the Georgi-Glashow model this 
is not the case, and the density of W is small. Thus, at Tatc there is a certain 
rearrangement of the thermal ensemble on short distance scales, but at long 
distances nothing happens — the string still confines. However, at Tc — Th/2 
the density of W reaches the critical value and the transition occurs. Note 
that at this temperature the large length fluctuations of the string are still 
suppressed — we are far below the Hagedorn temperature. The string is de- 
stroyed not due to "stringy" physics of the Hagedorn transition, but due to 
the short distance field-theoretical effects, viz. the fact that the fugacity of W 
is relatively large and that the interaction at short distances is Coulomb and 
not linear. 

9 The SU(iV) model 

So far in this review we have discussed the SU(2) gauge theory. Most of the 
discussion is easily generalized to the weakly interacting theory with SU(A^) 
gauge group, 

£ = -itrF^^i^'"' + trDf.^D^'^ - V{^) , (189) 

where 

$ = $"r°, (190) 

Here are traceless Hermitian generators of the SU(A^) algebra normalized 
such that tr(T"T'') = i<5°^ 

Depending on the form of the Higgs potential, there can be different pat- 
terns of gauge symmetry breaking. Since most of the details of the potential 
are unimportant for our purposes, we will not specify it except for restricting 
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it to the region of the parameter space where classically the gauge symnietry 
is broken to the maximal torus 

SU(iV)^U(l)^~i. (191) 

We also restrict ourselves to the weakly coupled regime, which means that the 
ratios Mw I are large for all N"^ — N massive W bosons. 

To characterize the perturbative spectrum of the theory it is convenient to 
use the Cartan-Weyl basis (i?\ E'^)i where W generate the Cartan subalgebra 
which is of dimension of the rank r = — 1 of SU(iV), 

[H\W]^0, [l,2,..iV-l], (192) 

and E" are the N{N — 1) ladder operators which satisfy 

[H\E^] ^a'E", (193) 

[E^, E^] ^N^0 E^+^ if a + P is a root (194) 

^2d-H if a = -/3. (195) 

The {N — l)-dimensional root vectors a — {a^ ,0^ , ...a^^^) form the dual 
Cartan subalgebra. There are obviously N{N — 1) such vectors corresponding 
to dim(SU(A^)) — rank(SU(A^)) but only A^— 1 of them are linearly independent. 
The non-vanishing inner products in the Cartan-Weyl basis are 

tT{H\ffi) = ^6'^, tiiE" , E^) = ^5"-^ . (196) 

At the classical level A^ — 1 gauge group generators are unbroken, which we 
choose to correspond to {H^}- Therefore, classically there are — 1 massless 
photons and A''(A'' — 1) charged massive W bosons. 

Our Weyl basis is chosen in such a way that the Higgs VEV is diagonal. 
Since the matrix <i> is traceless, there are A^ — 1 independent eigenvalues. In 
terms of the A^ — 1 dimensional vector h — {hi, ft,2, h^, ..h^^i) we have^ 

{^) = h-H, A^ = A^-H + Y,A^^E^. (197) 

a 

f For concreteness we order these numbers hi > /12 > .../ijv-2 > f'-N—iy which also breaks 
the discrete Weyl group. 
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For concreteness let us choose the following basis for the Cartan subalgebra: 
Hi = idiag(l, -1, 0, ...0), H2 = ^diag(l, 1, -2, 0...0) 

■■■ ^^-' = 71f^3ty'^'^^^^'^'^'-^'~^^~^^^- ^^^^^ 

As long a.s h ■ d for all roots, the gauge symmetry is maximally broken. 
The masses of the W bosons can be read off from the second term in the 
Lagrangian 

5'tr[^^, = y E A-'^AfhihjaW (199) 
=^Ms = g\h-a\. (200) 

The W bosons corresponding to the TV — 1 simple roots f3i, i = 1, N—1 
(an arbitrarily chosen set of linearly independent roots) can be thought of as 
fundamental, in the sense that the quantum numbers and the masses of all 
other W bosons are obtained as linear combinations of those of the fundamental 
W bosons. These charges and masses are 

= gp, = 9h-0. (201) 

As an example, consider the case of SU(3) broken down to U(l) x U{1). There 
are 6 massive W bosons. The simple roots can be taken as 

/3i = (^,^), /32 = (-^,^). (202) 

The remaining non-simple positive root is 

a3=/3i-/32 = (l,0). (203) 

The other three roots are — /Jj, —as. The masses of corresponding W bosons 

are 

Mwr = |(/ii + \/3/i2), Mw, = |(/ii - \/3/i2), Mw, = ghi (204) 
for hi > V3h2. 
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9.1 The monopole-instantons and the Polyakov effective Lagrangian 

Non-perturbatively, the most important contributions in the theory are due to 
the monopole-instantons. Those are classical, stable, finite action solutions of 
the Euclidean equations of motion arising due to the nontrivial nature of the 
second homotopy group of the vacuum manifold (n2(SU(A'')/U(l)^^^) = Z^). 
The magnetic field of such a monopole is long range. 



B. = ^9-H. (205) 

The A'^ — 1 dimensional vectors g are deteriiikifid by the non-Abelian general- 
ization of the Dirac quantization conditionE3Cj 

e'^'^-" = /. (206) 
Solutions of this quantization condition take the form 

N-l 

5 = -^n,/3* (207) 

where /3* are the dual roots defined by /3* — f3/\f3\'^. We will be working with 
roots normalized to unity, and thus fS* = j3. The integers rij are elements of 
the group 112 (see Ref. ^0|). The monopoles which have the smallest action 
correspond to roots taken once. The action of these monopoles in the BPS 
limit is 

Ms = —h ■ a = (208) 

9 9 

Just like with W bosons we can think of monopoles corresponding to simple 
roots as fundamental ones with magnetic charges and action 

47r ^ Att ^ -> 

gl = —A M, = —h-p,. (209) 

9 9 



For example, in the case of SU(3) (see Eqs. (|202[|20J)) the monopole action 
spectrum (in the BPS limit) is 

O'jr 27r 4-7r 

Mi = —{hi + V3h2), M2^—{hi-V3h2), Ah = —hi. (210) 
9 9 9 

The effect of these monopoles is to impart finite mass to all the perturba- 
tively massless "photons." The derivation of the effective Lagrangian follows 
exactly the same lines as the original derivation of PolyakovS for the SU(2) 
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theory. The resulting low-energy effective theory is wri 
N — 1 component field, fj, with the following Lagrangian: 



eff 



9 



327r 



, in terms of the 



(211) 



The sum is over all N{N — 1) non-vanishing roots. The potential induced by 
the monopoles is proportional to the monopole fugacity 



9' 



£,a = const. 



7/2 



9 



■ exp 



47rAfi 



9^ 



( Mh 



w 



(212) 



At weak coupling the photons are much lighter than the W bosons and thus 
are the only relevant degrees of freedom in the low-energy sector. 



9.2 The magnetic Zn symmetry 

The global symmetry structure is very important for the understanding of 
the deconfining transition. The relevant symmetry in the present model is 
the magnetic Zn symmetry. We now wish to explain how this symmetry is 
implemented in the effective low-energy Lagrangian — generalizing the result 
obtained for Z2 in the previous section. 

The order parameter of the magnetic symmetry is the set of magnetic vor- 
tex operators Vi, i — 1, ...,N — 1. These operators were constructed explicitly 
in Ref |l^. These operators carry the magnetic fluxes of the A'' — 1 U(l) 
Abelian magnetic fields. The defining commutation relation for Vi is 

— An 

[V,{x),Biy)] = mV,ix)5^{x - y). (213) 

9 

Here B is the A^ — 1 dimensional vector of magnetic fields]^ whose j-th compo- 
nent is the projection of the non- Abelian field strength onto the direction of the 
Cartan subalgebra generator Hj, and wj are A' — 1 weight vectors of SU(A^). 
The choice of the A^ — 1 out of N weight vectors is arbitrary. A change in this 
choice will lead to the redefinition of the vortex operators such that the new 
operators will be products of the old ones and their conjugates. It is always 
possible to choose these weights so that together with the "fundamental" roots 
Pi they satisfy the relation 

mPj = (214) 

''These magnetic fields can be constructed in an exnlicitlv gauge invariant way from the 
non- Abelian field strengths and the Higgs field, see Ref. |10|. 
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The flux eigenvalues in Eq. (213) are dictated by the requirement of the locality 
of the vortex operators and are analogous to the Dirac quantization co nditi on. 
The explicit form of the vortex operators in terms of the field 77 in Eq. (^ll|) is 



y.(f) = ^e'^' (215) 

VoTT 

with 

= • ?7 =^ rf = 2^/3-Xj. (216) 

i 

The effective Lagrangian can be written as a nonlinear a-model in terms of Vi 



KY.(v^v: - + ^ fc„ J] ^/-A (217) 

i a i 

with A — * 00. 

The matrix Aij = 2(3* ■ j3j depends on the choice of the fundamental roots. 
With the conventional choice of positive roots, where Pifij ~ —1/2, i ^ j, 
it is the Cartan matrix of the Lie algebra. All its diagonal elements are equal 
to 2, while all its off diagonal elements equal —1. We will find it however 
more convenient in the following to use a different set of fundamental roots, 
for which (3iPj = 1/2, i 7^ j. Such a choice is always possible for any SU(iV). 
For this choice of roots the off-diagonal matrix elements of Aij are all equal to 
1. 



For SU(3) we have 



and the effective Lagrangian 



A=( ^ 2 ) ^218) 



= d^Vid^V{ + ^V*d^ViV2d^V* + d^V2d^V* 

HiiViV^* + c.c) + ^2{V^V2 + c.c) + UViVi + c.c). (219) 

The magnetic Zn symmetry has an obvious and simple representation in 
this effective Lagrangian as 

Vi exp{27rm/7V}Vf . 
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As long as only small fluctuations of the phase fields Xi ^^'^ important, 
the Lagrangian ( pl9[ ) is equivalent to Eq. ( pll| ). Thus, at low temperature the 
descriptions based on these Lagrangians are equivalent. The difference appears 
just as for SU(2) only when the phase nature of Xi plays a role. Indeed, since 



Xi are treated in Eq. (219) as phases, dynamically one allows configurations in 



which these phases have nontrivial winding. On the other hand, in Eq. (211) 
such configurations cost an infinite amount of energy. As discussed in detail 
in Rcfs. |l^ and ^ the winding configurations correspond to the heavy W 
bosons. In fact, the explicit relation between the vorticity of the fields Vi and 
the electric charges is given byE3 

-mQ = / d/^5,,x;. (220) 
9 47r Jc^oo 

Thus, the difference between the two Lagrangians is important whenever the 
physics of the W bosons plays an important role. We have seen in the case 
of the SU(2) theory that W's are indeed important near the phase transition 
temperature. The same turns out to be true for arbitrary SU(iV). We iJiiis 
have to be careful to treat the W bosons properly in the transition regionEj 

9.3 The deconfining transition 

Close to the transition temperature we can safely use the dimensionally reduced 
Lagrangian. The zero Matsubara frequency sector is described by the two- 
dimensional Lagrangian 

= 32fe(^-^^)' + ^ 16^^"P^^" • (221) 

a 

However, as we noted before, our description should include W bosons, and 
so the fields ry should be treated here as phases with periodicity appropriate 



to Eq. (216). In fact, the Lagrangian also has to be augmented by a four 



derivative "Skyrme" term, whichJixes the energy of the winding states to be 
equal to the masses of W bosonsL3 Things are however simplified once we note 
that the density of W bosons at criticality is exponentially small due to the 
Boltzmann factor suppression. Thus, W^s can be treated in the dilute gas 
approximation in the same way as in the SU(2) theory. To do this explicitly 
we first have to understand how to wdte the partition function in the presence 
of one W boson of a particular typeEj 

Let us first consider a W boson corresponding to one of the fundamental 



roots /3fc. Combining Eq. ( |20l[ ) with ( |214[ ) and ( |220| ) we see that this W boson 



corresponds to unit vorticity of the field Vk and zero vorticity of all other fields 
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V}, j 7^ k. To create such a vortex in the path integral we must introduce an 
external "current" which forces the discontinuity of the field Xk 

Xk^Xk + 27r. (222) 

The partition function in the presence of one W boson is thus 

Z = J D[xix)] exp I - y d'y (^S[A] + ^ (a cos(2 ^ a ■ ftx.) j | , (223) 

where 

s[A] = Y^E^'^(^/^^' - Jliy,^))id^.x, - Jf.iy,x)), (224) 

where 

r^iy, x) = 27:d,kn^iy)S{y G C,) (225) 

with Cx a curve that starts at the location of the vortex (the point x), and 
goes to infinity, and is the unit normal to this curve. The insertion of this 
current forces the normal derivative of Xk to diverge on curve C, so that Xk 
jumps by 2tt across C. Since in the rest of the space Xk is smooth, the path 
integral is dominated by a configuration with unit vorticity of Xk ■ 



The path integral Eq. (224) differs from the partition function in the vac- 



uum sector by the linear term in the Lagrangian 

■J2 ^d^yP^■ Pj 3^X. = -/^ / dx^e^, djk ■ rf. (226) 



Defining in the standard way the dual field fj, 

id^,f] = €fj,^d^fl, (227) 

we can recast the contribution of this particular W boson in the form of the 
following extra term in the Lagrangian 

-*^/3l--^7- (228) 
This procedure can be repeated for a W boson corresponding to an arbitrary 



root a with the only difference that in Eq. (228), the root Pk is replaced by 
the root a. To create several W bosons one just inserts the external current 
which is the sum of the currents creating individual W's. 
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A dilute ensemble of such objects with small fugacities fia then has 

^ = nE^;^^^"/ l[d^^l[dy^zix,,yA (229) 



a n.ni 



The result of the summation over the number of W^s is the partition function 
with the Lagrangian 

^cff = g^^2y i^f^'^f + Caexp(m ■r])+Y^ M^exp {'^^:^^-'n] (230) 

a a ^ ^ 

with summation in both terms going over all non- vanishing roots of SU(A'^). 
The coefficients /Iq are proportional to the fugacities of the corresponding W 
bosons 

fia cx exp{-Mw^/T}. (231) 



Equation (23C) is the dimensionally reduced theory which we will now use to 
study the phase transition. 

Let us first disregard the W boson induced term in the effective La- 



grangian. If we do that, we are back to the theory of Eq. (221). An interesting 
property of this theory is that since the group is simply laced (all the roots 
are of unit length) the anomalous dimensions of all the interaction terms are 
equal. The scaling dimension of all the monopole induced terms is 

47rT , , 

Am = 232 
9 

This immediately tells us that at the temperature 

Tbkt = 1^, (233) 

all these interactions become irrelevant. Thus, at Tbkt one expects the Berezinsky- 
Kosterlitz-Thouless transition to take place. Above this temperature the in- 
frared behavior of the theory is that of iV — 1 free massless particles. Note 
that Tbkt does not depend on the number of colors N. If the picture just 
described were true, the universality class of the phase transition would be 
that of U^-i(l)- 

Of course, this is exactly the same situation as encountered in Ref. 25 in 
the SU(2) case. Again just like in SU(2) case this conclusion is incorrect due 
to the contribution of the W bosons. To see this, it is simplest to ask what 
would happen at high temperature if there were no monopole contributions at 



all. This amounts to studying Eq. ( pSOD with = 0. This theory describes 
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non-compact electrodynamics with iV — 1 photons and the spectrum of charged 



particles given by Eq. (201). This limit is again simple to understand, since 
the theory is exactly dual to the theory with monopoles and without charges. 
The scaling dimensions of all the W induced perturbations are equal and are 
given by 

Aw = (234) 

Thus, the perturbations are irrelevant at low temperature, but become relevant 
at 

Tjvc = |-. (235) 

Since Tjvc < Tbkt this tells us that we can not neglect the effects of charges 
at criticality. The story of SU(2) exactly repeats itself. Even the value of the 
temperature at which the scaling dimensions of the charge — and monopole 
induced perturbations are equal does not depend on N. 

We expect therefore that the actual transition temperature is 

Tc = ^, (236) 

at which point all perturbations have the same scaling dimension. This expec- 
tation is supported by the renormalization group analysis. 

The renormalization group equations for the theory Eq. ( |230| ) were studied 
in Ref. In general the equations are quite complicated due to the cross cor- 
relations between different operators. For this reason the space of parameters 
of the theory has to be enlarged if one wan ts to study the flow whose ultraviolet 
initial condition is provided by Eq. ( |23(]| ) with arbitrary values of fugacities. 
However, there is one simple case, that is when the initial condition is such 
that all the monopole fugacities are equal ^q.^ — — S,, and all the charge 
fugacities are equal = iiaj = This initial condition is stable under the 
renormalization group flow. In this subspace the renormalization group equa- 
tions, written in terms of the scaled temperature t = and dimensionless 
fugacities, read 

^=2n'Nt{^,^-e), (237) 
|^ = (2-i)M-2^(iV-2)M^ (238) 
1^ = (2-t)C-2^(iV-2)C'. (239) 
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These equations have exactly the property reflecting our previous discussion. 
That is the points t ^ 2, /i = and t = 1/2, £, = are both unstable. The 
stable infrared fixed point is 



^0 = 1 Mo = Co = — 7 r • (240) 

/ -,0 27r(Ar _ 2) ^ ' 

One can in fact easily check that, in the three-dimensional space of couplings 
t, ^ and fi^ this point has two attractive and one repulsive direction. This is 
precisely what one expects from the infrared fixed point located on the critical 
surface, the two attractive directions being the tangential directions to the 
surface. 

The renormalization group equations have an obvious duality symmetry, 
fj. S,, t ^ \/t. This is the reflection of the transformation r] fj on 



the level of the Lagrangian (230). The points i = 1, /i = ^ are symmetric 
under duality, and this ensures the existence of a self dual fixed point. This is 
important, since the exact position of the fixed point is scheme dependent. Its 
existence however is ensured by the duality symmetry. 

What is the nature of this fixed point? For = 2 we were able in Ref. ^ 
to fermionize the fixed point theory and show explicitly that it is equivalent to 
one massless Majorana fermion. We are not able to perform a similar analysis 
for arbitrary N. There are however several comments that we would like to 
make. Phase transitions in Zpf invariant spin models have been studied quite 
extensively. A recent discussion of the situation is given in Ref. ^ One 
considers a spin model of one phase field 9 with a symmetry breaking term of 
the type h cos{N9} which breaks the U(l) symmetry down to Z^. When the 
coefficient h of this symmetry breaking term is large, the model resembles the 
Potts model and thus (for > 4) has a first order phase transition. When 
the breaking is small on the other hand, the behavior is similar to the Villain 
model: the system undergoes two BKT type transitions with a massless U(l) 
symmetric phase at intermediate temperatures. At some particular "tricriti- 
cal" value of h the massless phase shrinks to a point and it comes together 
with the first order transition line. This tricritical point is self-dual and is 
described by a conformal invariant parafermionic theory with the central 
charge c = 2(A^ — 1)/{N + 2) introduced in Ref. In this type of model 
therefore one generically expects either a first order transition or a pair of BKT 
transitions with a massless phase in between. The tricritical behavior is special 
and requires fine tuning of the parameters. This is indeed also the prevailing 
general expectation for the order of the transition in 2-1-1 dimensional gauge 
theories at large N: either first order or Villain type U(l) invariant behavior. 

In fact, we find in our model a completely different situation. The tran- 
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sition is not first order, and there is no U(l) invariant massless phase. We 
stress that within the renormahzation group flow Eq. (239) the infrared fixed 



point Eq. (|24C|) has two attractive directions. This means that it governs the 
infrared behavior of the points which lie on a two-dimensional critical surface 
in the three-dimensional parameter space, and is therefore generic. This by 
itself does not preclude the possibility that this fixed point is the same as 
in the parafermionic Zj^ theory of Ref. If this is the case, it is quite 

interesting, since the point which appeared as "tricritical" from the point of 
view of usual spin models is in fact generic from the point of view of the 3D 
gauge theories. Although we can not prove that our critical point is described 
by the parafermionic theory, we will present some arguments supporting this 
conjecture. The point is that, as opposed to models considered in Ref. ^ 
our Lagrangian ( ^30| ) describes a theory of — 1 light fields. The theory of 
— 1 free massless fields has the ultraviolet central charge cuv = A^ — 1. 
However, this conformal field theory is deformed by the monopole and 14^- 
induced perturbations and flows to a different infrared fixed point. However, 
let us note that the central charge c = A' — 1 is precisely the central charge 
of the SU(Ar)i Wess-Zumino-Novikov-Witten (WZNW) model. The Ising (i.e. 
c = 1/2) model is the lowest among the minimal models with Virasoro (i.e. 
W2) symmetry. The highest model of this class is the c = 1 model (one free 
field ) which is precisely the SU(2)i WZNW model. When the c = 1 model 
was deformed by the monopole and VF-boson operators the central charge was 
reduced — and the resulting infrared theory was Ising-like. 

Now, Zn parafermions with c = 2(A^ — 1)/A^ -|- 2 are the lowest minimal 
models with Wn symmetry, and the highest is SU(A^)i (for more information 
about parafermions see, for example, Ref. ^ and references therein), which 
can be described in terms of A^ — 1 massless fields. Thus, if the theory in 
the ultraviolet describes A^ — 1 massless fields and has Wn symmetry, it is 
quite possible that result of the relevant (monopole-|-M^) deformation is a self- 
dual critical point. It is indeed known that the Z^ parafermion theory is the 
self-dual model with Wn symmetry. The fact that the central charge (and 
thus the effective number of degrees of freedom) is reduced in the process of 
the flow towards the infrared is of course in complete accord with Zamolod- 
chikov's C-theorem. It is therefore possible that the infrared fixed point that 
describes the universality class of the Georgi-Glashow model is the conformal 
Zn parafermion theory. 



Analysis of Ref. 54, although admittedly incomplete also supports the 
expectation that we do not have the Villain picture. As discussed in detail in 
Ref. |5^, it is the presence of the large number of fields that drives our theory 
away from the Villain behavior. 
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An interesting feature of this result is that the critical temperature in 
the SU(A^) theory at large N is proportional to the coupling and not to 
the 't Hooft coupling A = g^N. Thus, at large iV the critical temperature 
approaches zero. The physical reason for this is easy to understand. At large 
N and fixed A the Higgs VEV should also scale with N in such a way that 
the mass of W bosons remains fixed. The monopole action then grows as N 
and the photons get progressively lighter (exponentially with N)^. Thus, the 
thickness of the confining string grows and the density of W bosons needed to 
restore the symmetry becomes smaller and smaller. 

More importantly, our main conclusion is that the deconfining transition in 
the SU(A^) Georgi-Glashow model is second ord er and the universality class is 
determined by the infrared fixed point Eq. (240). This point is Zjv symmetric 
and self dual. We have given some arguments supporting the possibility that 
the fixed point theory is the Zn parafermionic modeO although we are not 
able to prove this explicitly. We can however definitely exclude the Potts and 
Villain universality classes. 



10 Concluding remarks 

We have discussed at length various aspects of the confining physics and the 
deconfining phase transition in weakly interacting gauge theories in 2+1 di- 
mensions. Naturally one should ask: where do we go from here? The ultimate 
goal we have in mind is of course strongly interacting 3+1 dimensional gauge 
theories. But this goal is still very far away. Many aspects of our discussion 
would change with dimensionality. In 3+1 dimensions monopoles are not in- 
stantons any longer, but particles, while the magnetic vortices are not particles 
but strings. The analog of the effective low-energy theories we have discussed, 
which incorporate the effects of these objects is not known. Moreover, weakly 
coupled phases of gauge theories are generically not confining and separated 
from the confining ones by phase transition. Thus, even if such effective the- 
ories were known it is not clear to what extent they can be used as a guide 
for the physics of the confining phase. The situation is not hopeless though, 
and a lot of work is being done on 3+1 dimensional theories. In particular, 
the dual superconductivity pictura^within the Abelian dominance hypothesis 
has been studied quite extensivelyEfl More recently the vortex-based ideas have 
also gained supportEl 

*This is analogous to the situation in QCD where the instantons become less relevant at 
large N and the rj' meson becomes massless. The major difference is of course that while 
the r]' mass in QCD decreases as 1/N, the photon masses in Georgi-Glashow model decrease 
exponentially. This difference is due to the non diluteness of the instanton gas in QCD as 
opposed to diluteness of the monopole gas in the Georgi-Glashow model. 
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Perhaps a simpler question is what can one say about the strongly inter- 
acting regime ir|_2+l dimensions. Here some pmgress has been made. Both 
monopole basedEJ and vortex based approachealj have been advanced. The 
basic encouraging feature is that the lightest degrees of freedom seem to be the 
same in the strongly and weakly coupled regimes.. This also seems to be the 
case near the deconfining transition temperatureLj Thus, a certain continuity 
is there and one hopes eventually to learn how to utilize it to the full extent. 

In short, the problems are hard — but don't despair ... 
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